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A CO J=3-2 map of M51 with HARP-B: Radial properties of the 
spiral structure 
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ABSTRACT 

We present the first complete CO J =3-2 map of the nearby grand-design spiral galaxy 
M51 (NGC 5194), at a spatial resolution of --600pc, obtained with the HARP-B instru- 
ment on the James Clerk Maxwell Telescope. The map covers the entire optical galaxy disk 
and out to the companion NGC 5195, with CO J=3-2 emission detected over an area of 
~9' X 6^(^21 X 14 kpc). We describe the CO J =3-2 integrated intensity map and combine our 
results with maps of CO J=2-l, CO J=l-0 and other data from the literature to investigate the 
variation of the molecular gas, atomic gas and poly cyclic aromatic hydrocarbon (PAH) prop- 
erties of M51 as a function of distance along the spiral structure on sub-kiloparsec scales. 
We find that for the CO J=3-2 and CO J=2-l transitions there is a clear difference between 
the variation of arm and inter- arm emission with galactocentric radius, with the inter- arm 
emission relatively constant with radius and the contrast between arm and inter-arm emission 
decreasing with radius. For CO J =1-0 line and Hi emission the variation with radius shows 
a similar trend for the arm and inter-arm regions, and the arm-inter-arm contrast appears 
relatively constant with radius. We investigate the variation of CO line ratios (J=3-2/2-l, 
J=2- 1/1-0 and J=3-2/l-0) as a function of distance along the spiral structure. Line ratios are 
consistent with the range of typical values for other nearby galaxies in the literature. The 
highest CO J=3-2/J=2-l line ratios are found in the central ~1 kiloparsec and in the spiral 
arms and the lowest line ratios in the inter- arm regions. We find no clear evidence of a trend 
with radius for the spiral arms but for the inter-arm regions there appears to be a trend for 
all CO line ratios to increase with radius. We find a strong relationship between the ratio 
of CO J =3-2 intensity to stellar continuum- subtracted 8 jam PAH surface brightness and the 
CO J =3-2 intensity that appears to vary with radius. 
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1 INTRODUCTION 

M51 (NGC 5194), the 'Whirlpool galaxy', is a nearby grand- 
design spiral galaxy viewed almost face-on (Table [T]), with two 
prominent spiral arms that can be traced close to the nuclear re- 
gion, and is undergoing an interaction with a small companion, 
NGC 5195, 4.5^ to the north-east. With its large angular size and 
strong spiral structure M 5 1 is well studied across a broad range of 
wavebands from X-ray to radio, and has been the subject of numer- 
ous studies investigating the origin of grand design spiral structure 



(e.g.|Tully||1974||Garcia-Burillo et al.|[T993bl [Shetty et al.||2007| 



[Meidt et al.||2008| ) and the relationship between spiral arms and 
star formation (e.g. Hitschfeld et al. 2009 , Foyle et al. 2010 ). It is a 
Seyfert 2 ( |Hoetal.|1997 ) and hosts a low-luminosity active galactic 
nucleus (AGN) surrounded by a ^100 pc disk of warm dense gas 
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( jKohno et a l. 1996 , Matsus hita et al.|1 998^. The two prominent spi- 
ral arms in M51, presumed to be caused via tidal interaction with 
the companion, are rich in molecular gas (e.g. |Aalto et aL]|1999[ 
[Shetty et al.|2007[|Koda et al.|2009| ), with molecular gas dominant 
over atomic gas until the outskirts of the disk (e.g. Garcia-Burillo| 
let al.|1993a||Nakai et al.| 1994] [Schuster et al.|2007| ). 

Molecular clouds provide the fuel for star formation, and in 
nearby galactic disks knowledge of the distribution of molecu- 
lar hydrogen gas is a crucial tool for investigating star formation 
and spiral structure. Since molecular hydrogen is difficult to detect 
directly, as it lacks a permanent dipole moment, carbon monox- 
ide (CO), observable at (sub)millimetre (submm) wavelengths, is 
widely used as a tracer of molecular hydrogen gas both in molecu- 
lar clouds i n the Milky W ay and in external galaxies. At a distance 
of 8.2 Mpc (^Tully|l9 74"l I" c orresponds to ~ 40 pc) and an inch- 
nation i = 20° (Tull y|l974] ), the proximity and face-on nature of 
M51 makes it possible to use single-dish (sub)mm telescopes to 
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study spiral arm structures and sub- structures at spatial resolutions 
of a few hundred parsecs. At this spatial resolution we are not sen- 
sitive to individual molecular clouds but rather clumps of molecular 
clouds, or Giant Molecular Associations (GMAs), as well as diffuse 
gas. The spiral arms of M51 contain numerous GMAs (as many as 
16 are noted by |Aalto et al.|1999| ) 

Since stars form not from the diffuse envelopes of molecu- 
lar clouds but from their dense molecular cores, tracing the warm 
dense gas is important for our understanding of star formation. The 
J =3-2 transition of CO (the J =3 level is 33 K above ground state) 
is a tracer of this warm, dense gas, while the lower critical density 
CO J=\-0 transition (the J=\ level is 5.5 K above ground state) is 
sensitive to the majority of gas in the cold ISM, including, for ex- 
ample, more diffuse gas. In fact, many authors have shown that the 
CO J =3-2 and CO J=l-0 lines cannot be tracing exactly the same 
components of molecular gas, and, importantly, that CO J=3-2 is 
the better tracer of warm dense gas (and hence star forming poten- 
tial) (e.g. [Wilson et al.|2009| and references therein). 

Although M51 has been mapped in low- J rotational ^^CO 
and ^^CO transitions in a number of single-dish and interferomet- 
ric studies, until now complete mapped observations of the full 
disk of M51 in the higher energy level transition CO J =3-2 line 
have been lacking, in particular because of the need for good tele- 
scope sites for higher frequency observations. Single-dish maps of 
the CO J=l-0 and CO J =2-1 transitions were presented by Garcia- 
Burillo et al. ( 1993a), Garcia-Burillo et al. ( 1993b), Schuster et al. 
( .2007; , .Nakai et al.. ( .1994; and .Koda et al.. ( .2009, .2011; (using 
the IRAM 30m and Nobeyama telescopes), and aperture synthe- 
sis maps have been presented by e.g. |Lo et al.|(1987|), [Sakamoto 



eTaT] ([1999}, [Aalto et al.| ([1999), [Regan et al.|p001| ) and JKodaet 



ar|([2009 2011). These studies have indicated that CO is mainly 
confined to the spiral arms. Other authors have studied tracers such 
as CO, CI, and CII in the central region or at individual positions 
in the spiral arms (e.g. [Israel & B aas "2002 1 [Kramer et al.[[2005l 
[Israel et aL][2006 ). CO J=3-2 observations to date have been re- 
stricted to the central regions and inner spiral arms of the galaxy 
(e.g. [Wielebinski et al.[ 19991 [Dumke ^t al.[2001| ) and have also had 
a relatively coarse angular resolution. Aperture synthesis mapping 
of the central region has also been carried out (e.g. Matsushita et[ 
[al.|2004[[Schinnerer et al._20 10), but aperture synthesis maps are 
known to resolve out and miss both a significant proportion of the 
CO flux and any extended emission. Single-dish CO J =3-2 obser- 
vations of the full disk of M51 are thus essential for determining 
the true distribution of molecular gas including any diffuse emis- 
sion in the inter-arm regions, and thus its relation to star formation 
and the total reservoir of star-forming gas. 

As well as understanding the distribution of molecular gas as a 
function of the spiral structure, it is also important to understand the 
relationship across multiple wavebands, from observations of the 
stellar component through to the dust and atomic gas components 
of the interstellar medium (ISM). Previous studies have suggested 
that polycyclic aromatic hydrocarbon (PAH) emission in the mid- 
infrared is a tracer of dust and gas in the cool ISM. For example, 
[Regan et"ar] |2006) found that within nearby galaxies PAH emis- 
sion is correlated with CO line emission, and several studies using 
Spitzer data have found that PAH emission is correlated with cool 
(~20 K) dust emission in nearby spiral galaxies (e.g. [Bendo et al.[ 
[2006l[Zhu et"aLl200 8 ). However, using HARP CO J=3-2 observa- 
tions of the spiral galaxy NGC 2403, [Bendo et al.[ ( [20T0l ) found the 
scale length of the CO J =3-2 radial profile to be statistically iden- 
tical to the scale length of the continuum- subtracted S^im (PAH) 
emission but found no correlation on sub-kpc scales. 



In this paper we present the first complete CO J =3-2 map 
of M51 covering the entire optical galaxy disk over an area of 
-- 10' X 7' (--24 X 17 kpc) at a spatial resolution of --600 pc. The pa- 
per is organised as follows. In Section [2] we describe our HARP- 
B CO J=3-2 observations, the data reduction and processing of 
the HARP-B data, and our treatment of ancillary data we took 
from the literature. Sections [3}|6] present our results. In Section [3] 
we present the HARP-B CO J =3-2 map, describe the distribu- 
tion of the CO J =3-2 emission in a qualitative way, and com- 
pare our CO J =3-2 map to maps of other tracers of molecular 
gas (CO J=2-l, CO J=l-0), atomic gas (Hi), star formation (Ha, 
24|^m) and PAH 8|am emission using data available in the litera- 
ture. We also investigate the contamination of submm continuum 
emission by in-band CO J=3-2 line emission. In Section[4]we dis- 
cuss how the appearance of M 5rs bright molecular gas rich spiral 
arms varies for different projected orientations and in Section [5] we 
disucss maps of the CO line ratio. In Section[6]we examine and dis- 
cuss the distribution and properties of the CO J=3-2 emission as a 
function of the distance along the spiral arms and inter- arm regions, 
and compare this to the distribution seen for other wave bands. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 HARP-B CO J=3-2 data 

Observations of the J=3-2 transition of ^^CO (rest frequency 
345.79 GHz) were carried out with the 16-pixel array receiver 
HARP-B ( [Buckle et aLl[2009] > at the James Clerk Maxwell Tele- 
scope (JCMT) in January 2007. The telescope beam size (Half 
Power Beam Width) at this frequency is ^\5" , which corresponds 
to ^600 pc at a distance of 8.2 Mpc. The 16 receivers of the HARP- 
B array are separated by approximately two beam widths {^30"), 
and a "basket-weave" scan pattern was used to achieve a final 
map sampling with spacing of half a beam-size. HARP-B was 
used in conjunction with the Auto-Correlation Spectral Imaging 
System (ACSIS) backend spectrometer, which was configured to 
have a bandwidth of 1 GHz and a resolution of 0.488 MHz (0.423 
km s~^at the observed frequency) per channel. 

The data reduction was carried out using Starlinljj software 
packages such as KAPPA, SMURF and CUPID in a similar man- 
ner to that described in detail in [Wilson et al.[ ( [2009| ) and [Warren[ 
[et al.[ j2010). Here, we give a brief overview of the data reduction 
process. First, the individual raw data files were inspected; obvi- 
ously noisy channels were removed, the data was despiked, and 
any bad baselines were flagged. Second, the raw scans were com- 
bined into a data cube using the task makecube that is part of the 
SMURF software package. In this process, a sinc(7r x)sinc(A:7r x) 
kernel was used to determine the contribution of of individual re- 
ceiver pixels to each map pixel (whose size was chosen to be 7.5^0; 
this "SincSinc" function is an often used weighting function that re- 
sults in a significant reduction in the noise in the map while having 
only a minor effect on the resolution. A baseline was then removed 
by fitting a fifth-order polynomial to line-free regions of the data 
cube. Thirdly, we created a mask of regions identified as having 
emission with a signal-to-noise >3 in a boxcar smoothed data cube 
(smoothed by 3x3 pixels spatially and over 25 velocity channels), 
using the isisk findclumps that is part of the CUPID software pack- 
age and uses the clumpfind algorithm ( [Williams et al.[1994j ). Lastly, 



^ The Starlink software package jCurrie et al.[[2008) is available from 
http://starlink.jach.hawaii.edu/[ 
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Table 1. Properties and observing parameters of M 5 1 



Parameter 




Value 


Centre R. A. (J2000) 




13^29^52.7^ 


Centre Decl. (J2000) 




+47°11M3'' 


Morphological type" 




SA(s)bc pec 


Distance ^ 




8.2 Mpc 


Position angle^ 




170° 


Inclination^ 




20° 


Vlsr'' 




+464kms-i 


Bandwidth 




IGHz 


Original channel width 




0.488 MHz (0.423 km s-i) 


rms(rX) per 0.423 km s 


■^ channel 


0.037 K 


rms (T^) per 10.6 km s~ 


^ velocity-binned channel 


0.008 K 


Beamsize 




15^^ 


Spatial scale 




600 pc/beam 








" de Vaucouleurs et al.|jl991) 




^ TuUy (e.g. 1974| 







we applied this mask to the original high resolution data cube and 
from this produced moment maps. The CO integrated intensity (ze- 
roth order moment) and velocity field (first order moment ) maps 
are shown in Figure [T] It is the CO J=3-2 integrated intensity map 
that is the focus of this paper. 

The size of the CO J=3-2 map is -10' x 7' (or 
^^24 kpcxl7 kpc), which covers the full galaxy disk. CO J=3-2 
emission is detected over the large proportion of this area (see 
Section [3]). The map rms noise level, measured over the main 
galaxy disk, is — 37mK (on the corrected antenna temperature, 
TJ, scale) for a channel with 0.423 km s~^ velocity resolution (or 
^-^8 mK for a 10.6 km s~^ channel). We note that pixels at the very 
edges of the map tend to exhibit an increased noise level, and we 
thus exclude these from our analysis. The la noise level in the 
total intensity map is estimated to be —0.1 Kkms~^(in TX units), 
based on the rms noise per 0.423 km s~^ channel and the number 
of channels at a typical position of emission. To convert the data 
from TX units to main beam brightness temperature (Tmb) units we 
divided by the main beam efficiency r]mb=0.6. The uncertainty on 
rjmb is 10-15 %. Throughout the paper, all results will be quoted in 
Tmb units. 

We compared our CO J =3-2 map to previous CO J =3-2 
Heinrich-Hertz Telescope (HHT) data: a map of the inner 3^-4' pre- 
sented by |Wielebinski et al.| ( |1999|) and a single HH T pointing at 
the central position presented by |Mauersberger et al.| ( |1999] ). After 
smoothing our map slightly to account for the difference in beam- 
size (the HHT beam-size was 22''and 18^' for the former and latter 
observations, respectively), we confirmed that both the shape of the 
spectra and the integrated intensity values in our map are in agree- 
ment with these previous studies within the calibration uncertain- 
ties (the uncertainty on rimh is 20% for the HHT data). We note that 
contrary to the findings of Wielebinski et al.| { [199 9), our HARP-B 
data does not show evidence of particularly extended CO J =3-2 
emission, and it appears that the HHT data may to some extent 
have been smeared out by an error beam. 

We also compared the integrated intensity values in the 
CO J=3-2 HARP-B map to previous JC MT CO J=3 - 2 obse rva- 
tions of the central regions presented by [Israel et aL| ( |2006| ) and 
find very good agreement. 



2.2 Ancillary data 

2.2.1 CO J=2-l and CO J=l-0 data 

For CO J=2-l and CO J=l-0, we used the integrated intensity 
maps presented in [Schuster et~ar] ( [2007| ) and [Koda et al.] ( [2009[ 
[201 1| ), respectively, which were kindly provided by the authors. 

The CO J=2-l dataset was observed with the HERA instru- 
ment on the IRAM 30-m telescope, whose beamsize at the fre- 
quency of the CO J=2-l line is r^\2" . The mean baseline rms is 
18 mK (TJ scale) at 5 km s~^ velocity resolution and the la noise 
level of the CO J=2-l map is 0.65 Kkms~\ on the TJ scale. To 
convert from TX units to main beam units we divide by r]mb=0. 57 
( [Schuster etal.[2007| ). For further details of the CO J=2-l observa- 
tions we refer to [Schuster et al.[p007[ ). 

The CO J=l-0 dataset was observed with the BEARS instru- 
ment on the Nobeyama Radio Observatory 45-m telescope, whose 
beamsize at 115 GHz is 15^^ The la noise level is 14.7 mK (TJ 
scale) in a lOkms"^ channel. To convert from TJ units to main 
beam units the data is divided by r]rnb=OA. For further details of 
the CO J=l-0 data we refer to [Koda et al.] ( [20091 [20Tl] ). 

As for the CO J=3-2 map, we also compared the CO J=2-l 
and CO J=l-0 maps to previous datasets and again find good 
agreement. 



2.2.2 Hi data 

We took H I integr ated intensity map s from The H I Nearby Galaxy 



We tooK n I integr ated intensity map s trom Ihe H I JNearby ualaxy 
Survey (THINGS ; [Walter et al.|2 008), via the THINGS websit£| 
THINGS was observed with the Very Large Array (VLA) in mul- 
tiple configurations to provide uniform sensitivity over a range of 
spatial scales, and provides two differently weighted datasets - a 
natural weighted integrated intensity map and a robust weighted 
map. In this work we use the natural weighted map as this should 
be more sensitive to extended, large-scale, emission and we do not 
require the higher resolution provided by the robust weighted map. 
The Full Width at Half Maximum (FWHM) of the Point Spread 
Function (PSF) of the H I image is 1 1.92'' x 10.01''. 



http : //www. mpia-hd . mpg . de/THINGS/O verview.html 
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2.2.3 Ha and Spitzer Space Telescope data 

We used H a (R-band subtracted), Spitzer Infrared Array Camera 
(IRAC) 3.6 and 8|j,m and Spitzer Multiband Imaging Photometer 
(MIPS) 24|j,m and 70|j,m data from the Spitzer Infrared Nearby 
Galaxies Survey (SINGS ;| Kennicutt et al.|2003^ ) sample, obtained 

For a description 



via the NASA/IPAC Infrared Science Archive 



of the data reduction technique we refer to [Regan et al.| ( |200 6). The 
FWHM of the PSFs are \.l" , 2", 6" and 18'' for the 3.6, 8, 24 and 
70|^m data, respectively. 

We used the 3.6 and 8^Lm images to produce a map of the PAH 
8|j,m emission using the relationship 

h{PAH 8|^m) = /.(8|^m) - 0.232/, (3.6|^m), (1) 



derived by |Helou et"ar] ( |2004| ), to subtract the stellar continuum 
from the 8|j,m image (where Ij^ has units of MJysr"^). For the 
two IRAC images used in this equation we first subtracted resid- 
ual backgrounds, by interpolating a smoothed version of the back- 
ground outside the optical disk, masked out bright foreground stars, 
and corrected to an 'infinite aperture' using values of 0.944 and 
0.737, respectively, as described in [Reach et al.| { |2005] ). 

For the 24|j,m image we also subtracted a residual background, 
using the same method of interpolating a smoothed version of the 
background outside the optical disk as we used for the IRAC im- 
ages. 



2.3 Alignment, convolution and re-gridding 

Since in this work we wish to make comparisons between physi- 
cal regions probed at different wavelengths, for example comparing 
structures in our CO J=3-2 image with CO J=2-l and CO J=l-0 
images taken from the literature (see Section [2.2.1] ), it is important 
to ensure the image coordinates are aligned and to convolve all the 
images to the lowest common resolution. We checked the relative 
pointing offsets between the three datasets by comparing peaks in 
the integrated intensity maps convolved to the HARP-B beamsize 
and pixel grid, by comparing the three datasets to the higher res- 
olution BIMA SONG CO J=l-0 map (Hei fer et al.|2003] ), and by 
examining the pixel-pixel relationship of the integrated intensity 
maps resampled onto a smaller pixel grid. We found no obvious 
offset between the two CO J=l-0 images, and also no obvious off- 
set between the CO J=l-0 images and the H a and 24|^m images. 
For the CO J=3-2 and CO J=2-l images we found offsets of +4!' 
and —A" , respectively, which is consistent with the pointing uncer- 
tainties for each dataset. A similar offset for the CO J=2-l dataset 
was previously reported by [Koda et al.[ ( [20T2] ). We applied these 
shifts to the HARP and HERA data accordingly. We stress that 
while small pointing offsets (a fraction of the beam) for each of 
the CO datasets is possible, it is unlikely that instead a single large 
(~8'') pointing offset could be present for just one of the datasets 
(and such a large offset in a single dataset would be highly obvi- 
ous). We note that although applying the pointing offsets results in 
generally better alignment of features in the two CO images, it does 
not significantly affect the overall results presented in this work. 

We then convolved all the images to match the PSF of the 
HARP-B data (15'' FWHM) and regridded them onto the same 
pixel grid as the HARP-B data (7.5" pixels). For the MIPS and 
IRAC images we obtained convolution kernels for the individual 
wavebands from [Gordon et al.]p008| ), which are produced using 



the PSF simulator STinyTim ( [Krist[[2002] ) and empirically deter- 
mined, respectively. Detailed information is given in [Gordon et al.[ 
( |2008[ ). From these, we created kernels that match the resolution of 
each dataset to the HARP-B CO J=3-2 PSF and appHed these to 
each dataset. For the CO and H I images we used Gaussian kernels. 
Ijin our analysis in Sections [5] and [6] only pixels detected at >3cr 
in the relevant images are used. 



2.4 The X-factor and the conversion of CO line intensity to 
molecular gas surface density 

The surface density of molecular gas, and related parameters such 
as the gas-to-dust ratio, are dependent on the ubiquitous CO-to- 
H2 conversion factor ('X-factor') Xco- Application of a constant 
X-factor across a galaxy assumes that there is little variation with 
metallicity or star formation activity. Yet, the metallicity is found 
to vary with radius in most spiral galaxies, leading many authors 
to suggest that Xqo may vary with metallicity in spiral galaxies 
(e.g. Arimoto et al.|1996 [ [Israel 2005 ), while other authors suggest 
that Xco may depend on other factors such as the cosmic-ray rate 
(e.g. [Bell et ^72006; . For example, [Garcia-Burillo erari ( [1993a[ ) 
found evidence that the value of Xqo may be a factor of about two 
smaller in the inter-arm region compared to the spiral arms. 

Moreover, the X-factor for conversion from CO line emission 
to molecular gas column density is usually based on the CO J=l-0 
line, and the value of Xqo appropriate for the J =3-2 transition 
of CO is rather uncertain. Previous authors have attempted to 
use CO J =3-2 data as an indirect measure of H2, by assuming a 
CO J=l-0 Xco factor and converting from CO J =3-2 intensity 
to CO J=l-0 intensity using the measured CO J=3-2/J=\-0 ra- 
tio. However, this method assumes that the Xco -factor is con- 
stant, and does not vary, for example, with galactocentric radius, 
and that the different CO lines trace the same molecular gas. Other 
authors have derived values for the CO J=3-2-to-H2 conversion 
factor, Xco(3-2)» using numerical models, for the central regions 
of galaxies where /co>100kms~^, and thus it is unclear whether 
it is applicable across the galaxy as a whole, in particular at much 
large galactocentric radii and in inter-arm regions where the CO 
intensities are considerably lower. 

Studies of the H2 surface density using the CO J=2-l and 
CO J=l-0 datasets used herein have been presented by previous 
authors (notably, [Schuster et al.[ ( [2007) and Nakai et al.[([l994)). In 
the present work, however, we are concerned with examining the 
variation of CO properties as a function of galactocentric radius 
and physical location (arm or inter-arm region etc.) across M51, 
and thus we choose to discuss the distribution in terms of observed 
CO integrated intensities and do not attempt to apply a CO J=3-2- 
to-H2 conversion factor. Also, examining relationships, e.g. the re- 
lationship between CO J=3-2 and PAH emission (see Section [6^ , 
in terms of CO intensities rather than molecular gas surface densi- 
ties allows us to make a consistent comparison with previous au- 
thors who disussed CO intensities (e.g. Regan et al.|2006| [Bendo[ 
[et al.[[2010| ). However, we will use the CO J=3-2 intensity to in- 
vestigate whether the CO J=3-2/J=l-0 ratio, or its dependence on 
another property of the galaxy such as CO J=3-2 intensity or galac- 
tocentric radius, could be used to scale the CO J=l-0 X-factor for 
use with the CO J=3-2 line. 



http://irsa.ipac.caltech.edu/data/SPITZER/SINGS/ 



Kernels are available from http://dirty.as.arizona.edu/~kgordon/mips/conv_psfs/conv_psfs.htm 
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Figure 1. Left: HARP-B CO J =3-2 integrated intensity (zeroth moment) map of M 51, on the T^b scale (Kkms~^), shown on a logarithmic colour scale. 
Pixels where no emission was detected are shown in black. The field of view (FOV) is ~10^ x7^ (or ~24 kpcx 17 kpc). The central position, listed in Table[T] 
is marked with a cross. Contours are at 0.5, 1, 2, 4, 7, 11, 18, 25, 32 and 39 Kkms~^; the \a noise level is 0.17Kkms~^). The peak integrated intensity 
is 41.3Kkms~^. The beam-size is ~15^^ and pixels are 7.5^^ in size (which corresponds to ~300pc at the distance of M51). Right: HARP-B CO J=3-2 
velocity field (first moment) map. Contour levels range from 380 to 580 kms~^ in intervals of 10 kms~^. 



A MAP OF M 51 IN CO J=3-2 EMISSION AND 
COMPARISON TO OTHER WAVEBANDS 



The HARP-B CO J=3-2 integrated intensity map of M51, shown 
in Figure [T] is the first sensitive CO map of the entire galaxy disk 
in the CO J=3-2 transition. The spatial scale is ^600 pc per 15^^ 
beam, or ^300 pc per pixel. The first contour level shown corre- 
sponds to 3cr, where a ^0.17 Kkms~^(Tmb)- CO J=3-2 emission 
is detected out to a radius of ~ 12 kpc, and is detected everywhere 
over a central area of ^14x 10 kpc. 



In this section we describe the CO J=3-2 image, which shows 
CO emission associated with molecular gas. We also give a qualita- 
tive comparison with images of other spiral arm tracers: lower CO 
transition tracers of molecular gas, H I emission that shows atomic 
gas, PAH emission (using the PAH S\im image produced as de- 
scribed in Section |2.2.3| l, H a that traces HII regions, and MIPS 
24|j,m emission that traces recent dust enshrouded star formation 
(and may also include a component of diffuse 24|^m emission that 
is not associated with recent star formation; Section [6|. These im- 
ages are shown at their native resolutions in Figures [TJ|3] with con- 
tours of the CO J=3-2 emission superimposed. 



3.1 The CO J=3-2 image 

The CO J =3-2 emission traces the two-armed spiral pattern, with 
one arm unwinding from the western side of the nuclear disk (here- 
after Arm I') and the other arm unwinding from the eastern side 
of the nuclear disk and stretching towards the companion galaxy, 
NGC 5 195, located at -- 1 1 kpc to the northeast of M 5 1 (hereafter 
Arm IF). The two arms are labelled in Figure [T] Arm I exhibits 
continuous CO J=3-2 emission out to a radius r ^250^^ (^10 kpc), 
with a notable 'kink' to the northwest at r ^110^^ where a peak 
of emission extends into the inter-arm region via a spur or bridge 
joining Arm I with Arm II. On the other hand. Arm II exhibits de- 
pression in the CO emission at ^100^^ west. This is not a void but 
rather the CO emission here is at a comparable level to that in the 
inter-arm region. 

Beyond this gap, there is a region of enhanced emission in 
Arm II at ^120^^ (^5 kpc) southwest, where the intensity is the 
highest of any emission in the spiral arm outside the central ~2 kpc. 
This emission peak is mirrored symmetrically opposite by such a 
peak in Arm I at 120^^ northeast. Enhanced star formation activity 
( ICalzetti et al.|2Q05] ) and [CII] emission ( Nikola et aL|2001| ) have 
previously been reported at these two locations, which are near the 
predicted corotation radius (160^^ Garcia-Burillo et al. T993a|b| l of 
the strong density wave pattern implied by kinematic studies (e.g. 
lAalto et aL|1999l|Garcia-Burillo et al.|19"93al|Kuno & Nakai|19971 
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Table 2. Typical values of CO J=3-2 line widths and brightness temperatures 



Central kpc Arms Outer arms Inter-arm 



Linewidth(kms-i) 40-80 [> 100] 20-50 15-20 15-20 

Brightness temperature (K) 0.3-0.7 [0.8] 0.1-0.5 0.1 <0.1 



|Shetty et al.|2007| ). Enhanced emission at these locations is inter- 
preted as being triggered by the interaction with the companion 
(e.g. [Toomre & Toomre„1972l [Nikola et ari[200T]l. As d escribed 
previously by a number of authors (e.g. |Nakai et al.|1994t , the lo- 
cations of these peaks also mark the locations beyond which the 
emission no longer follows a simple (e.g. logarithmic) spiral pat- 
tern. 

Along the inner part of Arm I the emission peaks are fairly 
regularly spaced (at intervals of typically ^1.8 kpc), while for 
Arm II fewer distinct peaks are discernible. The emission along the 
outer parts of both spiral arms appears to be more fragmented, as 
noted for CO J=2-l emission by [Schuster et al.| ( [2007| ). However, 
the outer part of Arm I to the west exhibits continuous low-level 
CO J=3-2 emission out to ^10 kpc (i.e. comparable to the level 
of inter-arm emission) with emission peaks fairly regularly spaced 
along the arm (at intervals typically of the order of ^2.5kpc). 

The arm-inter-arm intensity contrast is high, and the brightest 
CO J=3-2 emission clearly follows the spiral structure. Nonethe- 
less, inter-arm emission is detected at >3a everywhere out to a 
galactocentric radius of ^-^5 kpc and is detected for some regions 
out to ^10 kpc. The location of the highest arm/inter- arm contrast 
is at '^ 1.7 kpc south-west (just to the west of the brightest emission 
peaks in the image) where the intensity drops by a factor of ^14 
over just 9'\ We will return to the discussion of the spiral structure 
later in this section and in Section[6] 

There are several obvious spurs of CO J=3-2 emission con- 
necting the two arms. One of the brightest features is to the north- 
west, at a galactocentric radius of ^100^' (~4 kpc). To the south- 
west, at radii of ^50^^-110^' (~2-4 kpc), there are two further 
spurs, while another strong feature is seen to the north-east at a 
radius of ^150^^ (^6 kpc). These spur-like features are also very 
apparent in the CO J=2-l, CO J=l-0 and H I images, the latter fea- 
ture being especially apparent in the H I, H a, 24|am, and PAH S\im 
images (see Figures |2] and [3]). To the south, there is a hint of two 
faint spurs of emission connecting the two arms at ^-^3^-4^ (7- 
10 kpc). Although in these regions the CO emission is detected 
in typically only a handful of channels, with the signal-to-noise 
over these channels ranging from 2-7cr, we note that peaks in the 
stronger of the two spurs are also apparent in the CO J=2-l image 
(see Figure [2| and both features are clearly detected in recent Her- 
schel SPIRE 250|^m and 350|am observations ( Mentuch Cooper et| 
|al.|2012 ) which have comparable resolution. The CO J=3-2 image 
is seen to closely resemble the 250^Lm SPIRE image. 

Typical values of CO J=3-2 line widths and brightness tem- 
peratures at different locations in the galaxy are given in Table [2] 
As expected due to beam smearing, the broadest line spectra are 
found within the central kiloparsec (central ^-^25^^) of the galaxy, 
with typical widths of ^40-80 kms~\ and >100 kms~^ at the 
very centre. Here, brightness temperatures are found to be typi- 
cally 0.3-0.7 K, with the highest values (--0.8 K) found not right 
at the centre but corresponding to the locations where the spiral 
arms begin to wind out. The brightness temperature in the spiral 
arms is typically found to range from 0.1-0.5 K, with no obvi- 
ous dependence on radius until the outermost region of the arms. 



where values are ^--0.1 K. Typical line widths in the spiral arms, 
however, are found to decrease with distance along the spiral arms, 
from ^50 kms~^ to ^20 kms~^, as has previously been noted 
for CO J=2-l by |Garcia-Burillo et al.| ( [T993al ) and |Hitschfeld et al.| 
( |2009p who attribute the decrease in line width to streaming mo- 
tions (e.g. |Garcia-Burillo~ al. 1993b). In the outer spiral arms line 
widths are typically 15-20 km s~^, comparable to the widths found 
in the outer disk and in the inter-arm regions. Brightness temper- 
atures in the inter- arm regions are found to be <0.1 K. Overall, 
CO J =3-2 line widths in the centre of the galaxy and in the spiral 
arms are in very good agreement with those found for CO J=l-0 
by Na kai et al.|(|T994t and fo r CO J=2-l by [Garcia-Burillo et al.| 
( ,1993ar aM rHitschfeld et aT| ( [2009l ). 



3.2 Comparison to CO J=2-l and CO J=l-0 images 

CO J=2-l and CO J=l-0 integrated intensity images, at their native 
resolutions, are shown in Figure |2] All images are convolved to a 
common scale of ^600 pc and gridded on to the same 7.5^^ pixel 
grid as the HARP data, as described in Section [2] The CO J=2-l 
and CO J=l-0 emission is seen to trace the same two-armed spiral 
pattern as described for CO J =3-2. 

The CO J=2-l emission in particular closely traces the same 
features described above for CO J=3-2, with clear distinction be- 
tween the spiral arm and inter-arm regions. In fact, the peaks of 
CO J=3-2 and CO J=2-l emission are remarkably well aligned; 
at this spatial resolution we find no significant emission peaks in 
either map that are not in excellent agreement. Even in the cen- 
tral region where, at ^50^^ the spiral arms begin to wind out and 
Arm I is characterised by two peaks of CO J=3-2 emission, these 
two emission peaks are seen at the same spatial location in the 
CO J=2-l image. All the bright peaks of CO J=3-2 emission ap- 
pear to correspond to peaks in the CO J=2-l image and most of 
the bright CO J =3-2 peaks appear to have corresponding peaks in 
the CO J=l-0 image. Conversely, we find no bright peaks in the 
CO J=2-l and CO J=l-0 images that are not the also the locations 
of bright peaks in the CO J=3-2 image. 

The highest arm/inter-arm contrast in the CO J=2-l image is 
found at the same location as for the CO J=3-2 image, i.e. ^1.7 kpc 
south-west of the centre, although the CO J=2-l integrated inten- 
sity drops by only a f actor of about six over 9" , a nd falls by a 
factor of 14 over ^16^^ [Garcia-Burillo et al.||T993a| found a simi- 
lar result at this location using a different CO J=2-l dataset. These 
authors point out that the sudden drop between the spiral arm and 
inter-arm may in part be related to the fall-off of the nuclear com- 
ponent. 

The CO J=l-0 emission appears more diffuse with fewer ap- 
parent features. For example, the central CO J=l-0 emission is 
characterised by a single peak that is located between the peaks 
seen for the higher transitions. Indeed, a fundamental difference be- 
tween CO J=l-0 and the higher transitions is to be expected, since 
the CO J=l-0 line traces clouds of molecular gas, for example rela- 
tively low density cold diffuse gas, that are not (or only marginally) 
traced by the higher CO transitions. 
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Figure 2. CO J=2-l and CO J=l-0 integrated intensity images, from |Schuster et al.lpOO?} and |Koda et al.H2009||2011) , respectively, shown on logarithmic 
colour scales. Pixels where no emission was detected are shown in black. Images are shown at their native resolution (12^^and 16^^ respectively; see [2J but 
on the same 7.5^^ pixel grid as for the HARP image in Figure[T] Contours show the HARP CO J=3-2 integrated intensities as in Fig.[T] The FOV of all the 
images is identical to the HARP image shown in Figure[T] the spatial scale is indicated in Figure[T] 



For all three CO transitions the spiral arms show fairly regu- 
larly spaced concentrations or "clumps" of molecular gas, and, as 
noted for the lower CO transitions by previous authors (e.g. [Garcia^ 
[Burillo et al.|l 993 a[ [Schuster et al.,2007 ), at a spatial resolution of 
~600 pc may be referred to as GMAs and may represent either 
bound clusters or random superpositions of GMCs (e.g. |Rand &| 
|Kulkarni| 19901 IGarcia-Burillo et al.|1993a| ). The outer parts of the 
two arms appear more fragmented than the inner parts for all three 
CO transitions. 



3.3 Comparison to other wavebands 

The H a, 24|j,m, PAH 8|j,m and H I images are shown, at their native 
resolutions, in Figure [3] The PAH S^tm image clearly traces the 
clumpy dust structures of the ISM. The H I emission is stronger in 
the outer regions, and only a faint spiral structure is seen in the inner 
regions, as is typical in the inner parts of late type spirals where 
star formation is more active and the molecular phase dominates. 
By simple inspection of the images it is clear that the intensities of 
CO and H I are anti-correlated (also see Section |6!2] ). 

Compared to the CO J=3-2 emission, Ua, Hi, PAH 8|^m, 
and 24 \xm peaks in the northern parts of both spiral arms (within 
the co-rotation radius) are found to be slightly offset, by 300- 
400 pc (7"-\0''), towards the convex (downstream) side of the 
arms. This is particularly clear for Arm II at ^2.5 kpc north. This 



is in agreement with the findings of previous authors (e.g. 'Tilanus 
& Allen|199ri|Garcia-Burillo_eialJ1993a, Kunoerarir995 , Bas 



tian et al.||2005[|Egusa et al. 2009, Loui e, Koda & Egusa|[20T3] ), 
although we note that CO J =3 -2 (and likewise the CO J=2-l and 
CO J=l-0) emission is nonetheless present at the offset locations. 
Elsewhere in the spiral arms, even in the outer arms up to ^10 kpc, 
we see excellent agreement between bright peaks of CO J =3-2 
emission and bright peaks in the Hi, PAH 8|j,m, Ha and 24|xm 
images. One clear difference is that in the outer spiral arms of 
the CO J=3-2 (and likewise CO J=2-\) image, despite the good 
agreement between bright peaks we note a tendency for the Hi, 
PAH 8^Lm and 24|^m emission to be more extended on the convex 
side of the arms. This is particularly apparent at the western and 
north-eastern sides of Arm I, close to locations of spurs of emis- 
sion bridging the region between the inner spiral structure and the 
outer arms that are particularly clear in the PAH 8|j,m image. 



Comparison to the CO J=2-l image from [Schuster et al.| 
( [2007| presents a very similar result, since the CO J=3-2 and 
CO J=2-l peaks are coincident, and we find excellent agreement 
between bright peaks of CO J=2-l emission and all the other emis- 
sions even in the outer part of Arm I where these authors noted a 
slight offset. The relative locations of peaks in the CO J=2-l and 
H a images are in excellent agreement with those presented for dif- 
ferent CO J=2-l and H a datasets by IGarcia-Burillo et al.[([1993a|), 
and thus we believe that the slight offset noted by [Schuster et aL] 
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Figure 3. Images of (top row) Ho;, 24|Ltm, (bottom row) PAH 8|Lim (produced as described in Section[2j and Hi emission. Colour scales are chosen so as to 
best demonstrate the galaxy structure in each image (logarithmic colour scales for the PAH 8|Lim and 24 [xm images). For the H I image the colour scale shows 
the range from 0-0.3 Jy/beamkm s~^and, for the PAH 8|Lim image the colour scale ranges from 0-90 MJy sr~^ for M 51 (the companion galaxy NGC 5195 
has higher values). The images are shown at their native resolutions (Section|2] all have higher resolution than the CO images in Figures[^and |2j. Contours 
show the HARP CO J =3-2 integrated intensities as in Figure [T] 



A CO J =3 -2 map ofMSl with HARP-B: Radial properties of the spiral structure 9 




Position Angle (deg) 

Figure 4. Maps of CO J =3-2 (HARP-B, this work), CO J=2-l and CO J=l-0 integrated intensities shown in polar coordinates, where axes show galactocen- 
tric radius (in arcseconds; 1^^~40 pc) and position angle. The images have been de-projected using the parameters in Table^as described in Section [33] The 
images are shown on the same colour scales as in Figures [T] and |2] Contours for the CO J =3-2 image are the same as in Figure [T] For the CO J=l-0 image 
contours are at 2, 5, 10, 15, 20, 30, 40 and 50 Kkms~^ and for the CO J=2-l image contours are at 1, 3, 7, 18, 30, 41, 52, 64 Kkms~^(rmb). Arm I begins 
at RA.--190° and Arm II begins at RA.-^370° . 



( |2007| ) was likely due to a small pointing offset for which we have 
corrected in the present work (see Section|2]for details of our treat- 
ment of the datasets). 

The CO J =3-2 image also shows very good correspondence 
with with the PAH 8|^m image, with both images appearing to trace 
the same spiral spiral pattern. The locations of most peaks (and 
likewise troughs) in the CO J=3-2 image are also the locations of 
peaks and troughs in the PAH ^\im image. This is particularly ob- 
vious in the outer spiral arms and at the location of the spurs of 
emission spanning the inter-arm region to the southwest. The con- 



verse is also true, i.e. most features in the PAH 8|j,m image have 
counterparts in the CO J=3-2 image. There are, however, some 
clear differences. For example, to the south, where Arm I begins 
to unwind, there is a dip in the PAH 8|^m surface brightness that 
is not seen in the CO J =3-2 image. Also, in Arm II at ^5.5 kpc 
the CO J =3-2 peak is coincident with a peak in the PAH 8|j,m im- 
age, but there is a second PAH 8^Lm peak with no corresponding 
CO J=3-2 feature. The distribution of emission in the CO J=3-2 
and 24|j,m images is also generally similar, but the differences de- 
scribed for the PAH 8|^m image are also seen for the 24|j,m image. 
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Figure 5. As for Figure|4|but for H I, 24|Ltm and PAH 8|Ltm emission (from top). H I contours are at 36, 54,72, 108, 162, 234, and 288 mJy/beam km s~^ . 24|xm 
contours are at 1.5, 3, 5, 8, 15, 30, 65 MJy/sr. PAH 8M.m contours are at 0.5, 1.5, 3, 6, 12, 23 MJy/sr. 



One notable difference between the CO J =3-2 and 24|j,m images 
can be seen for the outer western part of Arm I, where the emis- 
sion appears very smooth at 24|j,m while several CO J =3-2 knots 
are seen along the arm. We will return the the subject of the rela- 
tionship between the CO J=3-2, PAH S\xm and 24 \im images in 
Section |6] 

Many of the same features can be seen across most wave 
bands, despite the different native resolutions. To allow a clearer vi- 
sualisaion of the comparison, and to aid the eye in delineating spiral 
features, in Figures [4] and [5] we display the CO, Hi, PAH 8|^m and 
24|j,m images in polar coordinates. These images have first been 
de-projected using the inclination angle, position angle and central 
position given in Table [T] 



In polar coordinates. Arm I can be seen to emerge from a ra- 
dius of ^50^^ in the central region at a position angle of ^190°. 
Arm II begins at -370° (i.e. 10°). Many of the features of the CO 
emission described in Section[3]are clearly emphasized in the polar 
images, and it is particularly striking how well-defined the spiral 
arm structure is in CO J=3-2 and CO J=2-l emission. Conversely, 
the CO J=l-0 image appears considerably more diffuse, with clear 
emission in the inter-arm regions. There are similarities, however. 
Spurs of emission connecting the two arms can clearly be seen in 
all three CO images at position angles of —60° and —340°, and 
the point at which the spiral pattern deviates from a simple loga- 
rithmic spiral is clearly discernible at position angles of —50° and 
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Figure 6. Percentage contamination at 850|xm by CO J=3-2 line emission, 
calculated as described in Section [34] Only pixels detected at >1.5cr in the 
CO J =3-2 and 850 |xm images are considered. Values are typically 10-20% 
in the spiral arms and are highest in the central ~2 kpc radius (on average 
~20% but reaching a maximum of 26% in the inner part of Arm II). In the 
inter- arm regions and outer spiral arms values are <10%. 



^-^230° . The latter are the locations of the symmetrically opposite 
bright peaks in each of the spiral arms described in Section [T7] 



3.4 Contamination of submillimetre continuum by CO J=3-2 
line emission 

Submm continuum observations of Galactic and extragalactic ob- 
jects may suffer from contamination from inband line emission, 
which for many sources is primarily CO J=3-2 line emission, 
due to the broad passband of many submm instruments such as 
SCUBA; the centre of the SCUBA 850|am bandpass is very close 
to the rest frequency of the CO J=3-2 transition (345.796 GHz). 
Continuum 850^Lm emission is typically used to trace cool dust, and 
since correcting 850|j,m emission for the contribution from CO will 
lead to a higher dust temperature and lower dust mass, knowledge 
of the amount of contamination can be vital. In particular, detailed 
knowledge of the amount of contamination by the CO J=3-2 line 
in a nearby galaxy such as M 51, where we can investigate the per- 
centage contamination as a function of radius and spiral structure, 
could be extremely useful for future galaxy studies. 

We carried out this investigation using our CO J=3-2 map to- 
gether with the total SCUBA 850|^m map from |Meijerink et al.| 
P005|), which includes emission from an exponential dust disk 



(see [Meijerink et al. 12005] for details and for a description of the 
data reduction procedure). Since both the CO J=3-2 and 850|j,m 
maps were obtained at the JCMT the beam sizes are compara- 
ble so we can compare the two maps directly and the only treat- 
ment of the 850|j,m map was to regrid it onto the same pixel grid 
as the CO J=3-2 map. Only pixels detected at >1.5cr in both the 
CO J=3-2 and 850^Lm images were considered. 

We first produced an image of the SCUBA equivalent flux pro- 
duced by the CO J=3-2 line at the recessional velocity of M51 
using the profile of the SCUBA filter bandpass and a 15^^ beam 
size. We then divided this image by the SCUBA 850|j,m image to 
produce an image of the percentage contamination at 850|j,m by 
inband CO J=3-2 emission. This is shown in Figure |6] Values are 
typically 10-20% in the spiral arms and are highest in the central 
^2 kpc radius (on average ^20% but reaching a maximum of 26% 
in the inner part of Arm II). In the inter- arm regions and outer spiral 
arms values are <10%. The average over the disk is 9%. We note, 
however, that if we perform the same analysis using the exponen- 
tial disk- subtracted 850|j,m map from [Meijerink et aL| ( |2005| ) the 
percentage contamination in the spiral arms is considerably higher. 



4 THE GEOMETRY OF CENTRAL MOLECULAR 
EMISSION 

Several authors have attempted to describe the geometry of the dis- 
tribution of molecular gas in the central regions of galaxies that are 
more highly inclined than M 5 1 based on the observed profile of 
the central gas component (often a 'two-horned' profile), leading 
to the conclusion that the gas is located in a molecular ring (e.g. 
[Irwin et al. 2011). However, the inclination of a galaxy may de- 
termine how reliably the geometry of its gas distribution can be 
determined. If, for example, a galaxy has bright molecular gas rich 
spiral arms, as is seen for M 51, this could potentially mimic a ring 
if observed edge-on (c.f. I srael et al. 2006 ). We can therefore use 
the M51 CO J=3-2, CO J=2-l and CO J=l-0 maps as a simple 
test of this scenario, by projecting the CO integrated intensity im- 
ages to several different edge-on orientations and examining the 
shape of the profiles along the major axis. 

Figure [T] shows major axis profiles of M51 for a major axis 
oriented along P. A. = 0, 45, 90 and 135°, produced using the pa- 
rameters given in Table[T] It can clearly be seen that the shape of the 
CO J =3-2 and CO J=2-l profiles at each projected orientation are 
similar. However, the important thing to note is that it is also clear 
that the shape of the edge-on profile can appear 'centrally-peaked' 
or 'two-horned' (see the second column of Figure |7] in particular) 
depending on the orientation. Thus it appears that, depending on 
viewing angle, a galaxy like M 5 1 with bright molecular gas rich 
inner spiral arms may be observed to exhibit a two-homed profile 
that mimics a molecular ring. 



5 CO LINE RATIOS 

Maps of the ratio of CO J=3-2/J=2-l CO J=2-l/l-0 and CO J=3- 
2/J=l-0 integrated intensities are shown in Figure [8] (hereafter 
i?3_2/2-i- ^2-1/1-0 and i?3_2/i-o)- These CO ratio maps were 
produced from the CO images shown in Figures [T] and [2] for pixels 
where the emission is detected above 3cr in both the input images. 
Input images (which are on the main beam temperature scale) were 
first matched to the HARP CO J=3-2 resolution and pixel grid as 
described in Section[2] 
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Figure 7. Profiles of CO integrated intensity (Tmb) along the major axis of M 51 when projected to various edge-on orientations (from left: major axis along 
P. A. = 0, 45, 90 andl35 degrees), using the parameters in Table[T] Rows show the different CO transitions as indicated. 



The CO ratio maps clearly exhibit a spiral-like structure, with 
the highest ratio values typically located in regions corresponding 
to the spiral arms and the lowest values typically located in the 
inter-arm regions. Values of the CO ratio in different regions of the 
galaxy are summarised in Table |3] 

We find a median R^_2/2-i ratio of 0.5 over the whole disk. 
In the central kiloparsec values are somewhat higher, ranging from 
0.6 to 0.7 at the central position and towards the location where 
Arm I begins to wind out. The ratio R^-2/2-1 clearly varies across 
the galaxy in a way that is similar to the spiral structure, with the 
highest values (^-^unity) found in the spiral arms and the lowest 
values (^0.1) found in the inter-arm regions. For i?3_2/i-o the 
median value is 0.36 over the galaxy disk and somewhat higher in 
the central kiloparsec. 

For the ratio of CO J=2-l/J=l-0 (i?2- 1/1-0) we find a median 
value of 0.80 over the galaxy and somewhat higher in the central 
kiloparsec. Similar results using the same datasets were presented 
As for i?3_2/2-i, the ratio R2- 



by 



Koda et al. 



(20121 



-1/1-0 



also 



appears to vary in a way that follows the spiral structure. However, 
the main difference is that, for i?2- 1/1-0 » the contrast between the 
spiral arms and inter-arm regions appears less pronounced over the 
central few kiloparsecs and the highest values are found in the spi- 
ral arms at larger galactocentric radii (>4kpc). While values of 
i?2- 1/1-0 at larger radii are found to be up to a factor of two higher 
than in the more inner spiral arms, values of i?3_2/2-i in the spiral 
arms appear to be similar at all radii. We note, however, that very 



high or low ratio values found in gradient locations such as at the 
edges of the spiral arms may be unreliable since, despite the care 
taken to align the maps and convolve them to a common beamsize 
(Section [23] ), there may still be minute alignment errors between 
the maps or minute differences in matched beamsizes. 

From inspection of Figure [8] while all three ratio maps appear 
to trace the spiral-arm structure, ratios with CO J=l-0 appear to 
show a clearer spiral pattern than i?3_ 2/2-1, which appears some- 
what smoother. This is not unexpected - since the higher CO transi- 
tions trace the spiral structure better than CO J=l-0, ratios of these 
transitions with CO J=l-0 will show a clear spiral pattern, while if 
the higher CO transitions trace the spiral structure roughly as well 
as each other then the ratio of the two will produce a map that is 
somewhat smoother in appearance. We discuss in more detail the 
variation of the CO ratios as a function of distance along the spiral 
arms and inter-arm regions in Section[6] 

We find the values of R^_2/2-i to be in agreement with the 
range of values found for M51 by previous authors and for spi- 
ral galaxies in general (e.g. [Mauersberger et alT]|1999| ). The aver- 
age i?3_2/i-o line ratio over the whole disk of M51(~0.36) falls 
within the range of ratios found by 'Mauersberger et al.| { [1999| 0.2- 
0.7), and is consistent with the average values found in a number of 
other recent CO J=3-2 studies of nearby galaxies with the HARP- 
B instrument (e.g. [Wilson et al.|2009| [Warren et aL,2010p . These 
authors point out that if galaxies have enhanced line ratios in their 
centres, then previous observations such as those by [Mauersbergerj 
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Figure 8. Maps of the ratio of the J=3-2/J=2-l, J=2-l/J=l-0, and J=3-2/J=l-0 CO integrated intensities (from left), shown on colour scales that are centred 
on the median ratio value in each map and have a range of la from the median. Median ratio values are 0.5, 0.80 and 0.36, respectively (see Table [3j. The 
maps are produced as described in Section[5] the resolution, pixel size and FOV of each input image has been matched to those of the HARP CO J=3-2 image 
(Figure [^, and only pixels with emission >3cr in both input images are shown. The uncertainties in the CO ratios are mainly dependent on the calibration 
uncertainties for the individual CO images. 



ITable 3. CO line ratios in M 511 



Global 

(1) 



Central kpc Arm Inter- arm 

(2) (3) (4) 



COJ=3-2/J=2-l 
COJ=2-1/J=1-0 
COJ=3-2/J=l-0 



0.50t ±0.14* [0.1-1.0] 

0.80§ ±0.32* [0.3-3.1] 

0.36* ±0.17* [<0.1-1.7] 



0.63 


0.5 


0.4 


0.86 


0.9 


0.5 


0.54 


0.4 


0.2 



(l)-(2) Median CD ratio over tiie tuU galaxy disk and central kiloparsec, tor pixels where the emission is detected above 3cr 
in both input images; the uncertainty indicates the standard deviation of the ratio over the whole disk. 

(3)-(4) Typical CO ratios in the spiral arms and inter-arm regions. 

fThe median "CO J =3-2 line intensity is 3.6 K lkms~^h nd the median ICO J=2-ll line intensity is 7.5 K'kms~M 
*The standard deviation of the CO line intensities in the input images is 6, 10 and 5 for ICO J=3^2.iCO J=2-T1 and lCO J=1-0I 
respectively. The sources of uncertainty on the CO line intensities are described in Section[2] 
§The median ICO J=2-ll and lCO J= 1-0 line intensities are 7.6 and 5.5 K^^^^ respectively. 
[|]rhe median lcb J=3-2l and lCO J=Fol line intensities are 2.2 and 3.7 K lkms~^l respectively. 



|et al.|(T999] ) that only observed the central regions of galaxies may 
have overestimated the ratios averaged over the whole galaxy disk. 
The factor of >10 variation in Rs-2/1-0 across the disk of ]V[51 
is high compared to some previous observations (e.g. |Thornley &| 
|Wilson[1 994, JVLauersberger et al. 1999) but is in good agreement 
with the variation found in recent studies by Bendo et al.| ( |2010| l 
and [Wilson et al.| p009| ) for other nearby galaxies observed in 
CO J=3-2 with HARP-B. As these authors suggest, a large vari- 
ation in i^3_2/i-o may be real but was not observed before due to 
some regions being undetected or unresolved. 

In the present work, no attempt has been made to extract 
molecular gas parameters from the measured line ratios, because 
the observed ^^CO transitions are strongly degenerate with respect 
to molecular gas temperature and density. This degeneracy can be 
broken by e.g. observations of an optically thin isotope such as 
^^CO, but these are not available at present on a galaxy- wide scale. 



Relevant discussions of a few individual pointings may be found in 
[Kramer et all ( [2005] ) and jlsrael et aL| ( [2056l ). 



The nuclear region of IM 51 appears to have line ratios that are 
different from those of much of the disk, indicating perhaps an ex- 
citation mechanism substantially different from the presumed PDR 
excitation of the disk gas (also see I srael et al.^2006j . Whether or 
not this related to the low-luminosity AGN referred to in the intro- 
duction is not presently clear. To answer this question reqwuires a 
more detailed study of the molecular line emission from the central 
region. 
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6 RADIAL PROPERTIES OF M 51 AS A FUNCTION OF 
THE CO SPIRAL STRUCTURE 

In this section we investigate the properties of M 51 as a function of 
galactocentric radius and, specifically, as a function of the distance 
along the observed molecular spiral structure. 



6.1 Describing the CO spiral structure 

The commonly used method of obtaining radial profiles by averag- 
ing in elliptical annuli over the whole disk is likely to average out 
the true variation of properties along spiral arms and in inter-arm 
regions, so here we describe our alternative approach of producing 
separate radial profiles of the spiral arms and inter-arm regions. 

The spiral pattern at small galactocentric radii (R< 120^0 in 
M 5 1 has often been described with simple symmetrical logarith- 
mic spirals (see e.g. Nakai et al.__1994, Shetty et al.|2007| ). How- 
ever, the spiral arm pattern clearly no longer follows a logarithmic 
spiral at R>120^^ (>5 kpc), with the fracture in the spiral pattern 
at R~120'' and the distortion at radii beyond this presumably due 
to the tidal interaction with the companion galaxy; in the images 
in Figures |4] and [5] this can be easily seen at RA.^50° for Arm I 
and RA.'-^230° for Arm II, for all three CO transitions as well as 
PAH 8^Lm and 24|j,m emission. In previous studies tracing the spi- 
ral structure in M51 to radii much larger than 120^^ authors have 
adopted a variety of approaches for defining the spiral structure, but 
have often involved selecting regions based on CO surface bright- 
ness. For example, |Hitschf eld et al.| ( [2009| ) attributed all emission 
above a certain threshold molecular gas surface density to spiral 
arms and the remainder to the inner-arm regions, and |Foyle et al.| 
P010| ) followed a similar approach but instead used the surface 
density of stellar 3.6|am emission. This approach, however, has the 
effect of including emission features (spurs) bridging the gap be- 
tween the spiral arms as part of the arms themselves. 

In the present work, we instead wish to examine the proper- 
ties M51 as a function of physical location and distance along the 
observed spiral morphology out to where the interaction with the 
companion NGC 5195 is taking place (^11 kpc), and thus we re- 
quire a different approach. We take the simple approach of defining 
a mask by eye that broadly describes the observed appearance of 
the CO spiral structure out to large radii - i.e. two spiral arms and 
two inter-arm regions, where we define the "spiral arm" region as 
a continuous broad structure that is linked seamlessly from the nu- 
clear region in a broadly spiral pattern. The remaining region in 
between the spiral arms we refer to as the "inter-arm" region, ir- 
respective of any emission features it may contain. The resulting 
mask - defined based on the CO J=3-2 image while ensuring that 
it adequately describes the spiral structure observed in the other 
two CO transitions - is shown in Figure [9] The width of the spiral 
arms defined in this way is ^1 kpc (^20^^-30^0, when corrected 
for the telescope beam, which is in excellent agreement with the 
values found for CO J=2-l and CO J=l-0 by |Garcia-Burillo et al.| 
( |1993a| l and |Nakai et"aL] ( |1994 |), respectively. 

We produce radial profiles by plotting the pixels in each spiral 
arm or inter-arm as a function of galactocentric distance (where the 
central position is defined as given in Table [TJ. These are shown by 
the black points in Figures [To| (filled and open points show Arms I 
and II, respectively) and|ll| We also bin this data into 0.5 kpc 
bins (approximately the size of the HARP-B beam- size), and this 
is shown as the red (Arm I) and blue (Arm II) points in the figures. 
The latter is effectively the equivalent of averaging the emission in 
a series of circular annuli separated by 0.5 kpc in which only J^- 
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Figure 9. Mask used to describe the spiral structure (i.e. used to differentiate 
arm and inter-arm regions) in M 51 (shown in green), derived from the CO 
(J=3-2 and J=2-l) datasets as described in Section [6| The greyscale image 
shows an H a image projected to a face-on orientation using the parameters 
in Tabled 



tected pixels are considered. We do not set undetected pixels to zero 
since there is considerable difference in the sensitivity of the three 
CO maps we wish to compare (see Section|2]). We note that the CO 
profiles can be reliably interpreted out to ^5-6 kpc, beyond which 
the number of detected pixels is relatively small. 

All radial profiles were produced using de-projected images 
produced as in Section |3.3| and that emission within the central 
1.7 kpc is not considered for the analysis of the spiral structure. In 
the case of the CO ratios, we only consider pixels that are detected 
above 3 a in both input CO maps. 



6.2 Radial variation of CO, H i and PAH emission 

Figure [Tol shows the radial profiles of the CO J=3-2, CO J=2-l, 
CO J=l-0 and H I integrated intensities. The CO profiles of the two 
spiral arms are generally similar, with CO intensity decreasing with 
radius but with a strong 'humps' at ^5 kpc and ^8 kpc. The Hi 
profiles also show some evidence of a broad humps at these radii. 
The former is the radius at which two strong regions of CO emis- 
sion are seen in the spiral arms, to the north-east and the south-west 
of the centre. Two notable differences are that the scatter appears 
to increase somewhat towards higher CO transitions, and there ap- 
pears to be less scatter for the Arm II than Arm I. This is probably 
related to the interaction with the companion to the north of M 5 1 . 
Conversely, there is remarkably little scatter for the inter-arm emis- 
sion, with the shape of the profile broadly following that seen in 
the spiral arms (but without the broad hump). This is what may 
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Figure 10. Radial profiles of velocity-integrated CO J=3-2, J=2-l, J=l-0 and H i intensity and PAH 8|xm surface brighness in the two spiral arm and inter- 
arm regions (in units of K km s~ ^ (T^]^) and Jy km s~ ^ , respectively). Arm and inter-arm regions are defined using a mask based on the CO J=3-2 emission 
that is modified so that it also adequately describes the spiral structure observed in the other wavebands (see Section[6|for further details); emission within the 
central 1.7 kpc is not considered. Black points show individual 7.5^^pixels with a signal-to-noise >3; filled points indicate Arm I and open points indicate Arm 
II (see text for definition). Red and blue points show the same data binned into 0.5 kpc bins (approximately the size of the HARP-B beam-size), for Arm I and 
Arm II, respectively. Error bars show the error on the mean, and where not present are smaller than the symbol size. The third column shows the arm-inter-arm 
contrast as a function of radius for the two spiral/inter-arm regions, i.e. the ratio of the first two columns; see Section[63]for definition). 
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Table 4. Properties of the exponential disk in CO J=3-2 and CO J=2-l compared to the properties of the total gas exponential 
disk presented by Hitschfeld et al. (2009 1 and the exponential dust disk presented by ,Meijerink et al., ( ,2005j . 



Exponential disk determined in: 


Scale length 


disk/tota 


mass fraction 




(kpc) 




% 


Molecular gas (CO J=3-2) 


6.9±0.4 




14 


Molecular gas (CO J=2-l) 


5.3±0.1 




25 


Total gas (CO J=2- land Hi) 


7.4 




61 


Dust 


6.2 




55 



be expected if the inter-arm emission is diffuse emission underly- 
ing the disk. Despite these differences, the total profile (Figure [T0| 
for the whole galaxy disk shows a strong similarity between the 
shapes of the binned CO profiles on 0.5 kpc scales. The emission 
scale lengths we measure for CO J=3-2 and CO J=2-l are shown 
in Table [4] For comparison we also show the scale lengths of the 

1 ([20091 ) and 



Table 5. Average arm-interarm contrast. 



total gas exponential disk presented by [Hitschfeld et al 



the exponential dust disk presented by'Meijerink et al.|i 



20051 



The PAH8|j,m surface density profile is seen to be similar 
to the profiles of CO J=3-2 intensity but is very different to the 
profile of Hi emission, suggesting that in M51 PAH emission is 
not associated with H I. While the radial profiles of CO J=3-2 and 
PAH 8|am emission are similar to one another, the radial profile 
of the Hi emission is clearly very different (Figure [TOl, as is also 
obvious from simple inspection of the images in Figure |3] As pre- 
viously noted in Section |3.3| it is also clear from Figure [To| that 
the intensities of CO and H I are anti-correlated, i.e. that the two 
species are complementary, which suggests that the ratio H 1/H2 is 
a strong, and apparently continuous, function of radius. 



6.3 Radial variation of arm-inter-arm contrast 

In the right-hand panel of Figure [T0| we show the contrast between 
arm and inter-arm emission as a function of radius for the two spi- 
ral arms and inter-arm regions. Contrast is defined as the ratio of 
the mean spiral arm velocity-integrated intensity versus the mean 
inter-arm integrated intensity in bins of 0.5 kpc (i.e. the ratio of the 
red points shown in Figure[To]). The average arm/inter- arm contrast 
values for CO J=3-2, CO J=2-l, CO J=l-0 and Hi emission, as 
well as PAH 8|^m surface density, are given in Table|5] As expected 
from simple inspection of the images in Figures [TJ|2[there appears 
to be a trend for arm-interarm contrast to increase towards higher 
CO transitions. For CO J=2-l, a value of 2.9 was previously found 
by [Hitschfeld et al.] ( |200 9) who used used the same dataset as the 
present work but with different method of defining the arm and 
inter-arm regions, while we find a value of 2.1. Both values are in 
reasonable agreement with |Garcia-Burillo et al.| ( [T993b| ). 

It is clear that for Hi and CO J=l-0 emission there is little 
variation in the arm-inter- arm contrast with radius over ^-^10 kpc. 
Conversely, for the higher CO transitions the contrast is high in the 
central few kiloparsecs and clearly decreases towards larger radii. 
The trend is particularly clear for CO J=3-2 emission, which de- 
creases from a contrast of '^ 4 at 2 kpc to ^-^1 at ^6.5 kpc (the 
co-rotation radius). 



6.4 Radial variation of CO line ratios 

we show the CO ratios Rs-2/2-1^ ^3-2/i-o» and 



Transition 


Arm-interarm contrast 


CO J=3-2 


2.4±1.1 


CO J=2-l 


2.1±0.5 


CO J=l-0 


1.7±0.5 


Hi 


1.6±0.3 


PAH 8|Lim 


2.2±0.8 



11 



In Figure 

i?2 -1/1-0 as a function of distance along the spiral arms and inter- 
arm regions. The last column shows the radial profile for the entire 



galaxy. Despite considerable scatter in the individual pixel values, 
averaged on 0.5 kpc scales the variation in CO ratio in the spiral 
arms is surprisingly small, typically much less than a factor of two, 
both within an individual spiral arm and compared between the two 
spiral arms. The main notable difference between Arm I and Arm II 
occurs for Rs-2/1-0 at ^7 kpc, near the co-rotation radius. Here, 
the value of i?3_2/i-o for Arm I decreases by a factor of about two 
while for Arm II the ratio increases by a similar factor. While the 
values for the other CO ratios also increase for Arm II at this radius, 
for Arm I virtually no significant CO J=2-l emission is detected 
beyond ~6 kpc. The CO ratio profiles for the total galaxy disk again 
emphasize that the variation on 0.5 kpc scales appears to be within 
a factor of two, whether in individual spiral structures or averaged 
over the whole galaxy. It is interesting to note, however, that while 
the total profile is almost flat from 2-4 kpc in all three transitions, 
in the individual spiral arms there is nonetheless some variation - 
in the case of i?3_ 2/2-1 as much as a factor of two. Comparing 
Arm I and Arm II it is clear why this is - in the eastern arm the 
ratio is decreasing while in the western arm it is increasing. Thus, 
this shows that in only taking radial averages over the whole galaxy 
one would 'miss' a factor of two variation between the two spiral 
arms. 

Median line ratios are i?3_2/2- 1^^0.5 in the spiral arms and 
^-^0.4 in the inter-arm regions. For i?2- 1/1-0^ average values are 
^^0.4 in the spiral arms and ^^0.3 in the inter-arm regions. We note 
that we find lower median values of i?2- 1/1-0 and i?3_ 2/2-1 for 
inter-arm I than for inter-arm II. While this may hint at a small 
systematic difference in the average CO ratio between the arm and 
inter- arm regions we note that the difference is only of the order of 
one standard deviation (Section |5]). Lower line ratios suggests that 
the molecular gas may on average be cooler and/or less dense, so 
it is not unexpected to find higher values in the spiral arms which 
are clearly the sites of star formation (see e.g. Figures [4] and [5]) 
and thus should host warmer, denser, gas. Our values of Rs-2/2-1 
and i^2- 1/1-0 are somewhat lower than found in previous works - 
[Garcia-Burillo et aL] (|1993a) and Kramer et al. (2005) find values 
of ~0.7 over the whole disk and in the central region, respectively. 

However, as can be seen in the last panel of Figure [TT] our 
value of i?3_2/2-i in the central kiloparsec is ^-^0.6, with values 
for individual lines of sight in the inner spiral arms ranging up to 
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Figure 11. CO line ratios as a function of distance along the spiral arms and along the inter-arm regions, starting from a radius of 1.7 kpc. Symbols are the 
same as for Figure[To| Median line ratios are listed in TablejS] The final column shows the total radial profile for the galaxy, including the central 1.7 kpc. CO 
emission is detected everywhere, in both the relevant transitions, out to a radius of ~4kpc. Thus, out to ~4kpc the binned profiles shown here are equivalent 
to the profile that would be obtained by averaging in a series of concentric annuli. 



^0.8, while our value for i?2-i/i-o are '-^O.S-O.G in the central 
kiloparsec. Thus, within the uncertainties, we find our values of 
i?3_2/2-i to be in agreement with those in the literature. Our val- 
ues of i^3_2/i-o are also in good agreement with those found in 
previous works. We find median values of ^0.4 in the spiral arms, 
~0.2 in the inter-arms, and ^0.54 in the central kiloparsec. For 
comparison, [Garcia-Burillo et al.| ( |1993al ) found values of ^0.6 in 
the spiral arms within the central l^ 

For each of the CO line ratios, we generally find no significant 
difference between the mean values in Arms I and II, nor do we find 
any significant difference between the two inter-arm regions. Thus, 



while there is clearly variation between the individual spiral arms 
at a given radius, as described above, when averaged over distance 
the CO ratios appear to be similar. 



The ratio i?2-i/i-o for the spiral arms shows some evidence 
of an increasing trend towards larger galactocentric radii, while 
there is no clear trend for i?3_ 2/2-1 or i?3_ 2/1-0- The lack of a 
clear trend for i?3_2/i-o to vary with radius when studied over ei- 
ther all regions of emission, (the right-hand panel of Figure [TTJ or 
just the spiral arms (the left-hand panel o f Figure pT| ), is in agree- 
ment with the findings for NGC 2403 by |Bendo etal.| ( |2010| ). How- 
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Figure 12. CO J=3-2/J=2-l, J=2-l/J=l-0, and J=3-2/J=l-0 line intensity ratios (rows) plotted as a function of CO J=3-2, CO J=2-l, and CO J=3-2 line 
intensity, respectively, for the spiral arms (left panels), inter-arm regions (centre panels) and entire galaxy disk (right-hand panels). Only pixels with emission 
detected above 3cr in both CO transitions are shown. In the right-hand column we colour-code the plotted points according to galactocentric radius as indicated 
in the top right panel: < 1 .7 kpc (red), 1 .7-3 kpc (green), 3-4.5 kpc (cyan), 4.5-6 kpc (blue), and >6 kpc (purple). For the left and centre panels, filled and open 
points show Arm I and II, respectively. 



ever, for M 51, in the inter- arm regions there appears to be a trend 
for all the CO ratios to increase with radius. 



6.5 CO line ratios and the X factor 

The large pixel-to-pixel variation in the CO ratio i?3_ 2/1-0 de- 
scribed above, together with the suggested different behaviour with 
radius for the spiral arms and inter-arm regions, implies that ob- 
taining molecular gas column densities from the CO J =3-2 data 
by applying a simple scaling term to match it to the X -factor de- 



termined for CO J=l-0 is likely to introduce large uncertainties 
for sub-kiloparsec scales. As in |Bendo et al.| ( [2010l l, we investigate 
this further by looking at the relationship between CO line ratio 
and CO surface brightness. We plot this for Rs-2/2-1^ ^2-1/1-0 
and i?3_2/i-o in Figure 12 In the first two columns we plot sepa- 



rately the regions of spiral arm and inter-arm emission, as described 
for the previous figures, and in the last column we show the entire 
galaxy disk (where emission is detected at >3a in both CO lines), 
colour coded for galactocentric radius. We can immediately note 
several features of the distribution of the CO ratio with CO surface 
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Table 6. Parameters of linear fits in Figures[T3]and[T8] 

Radius A B Vp 

(kpc) 

Relationship between CO J=3-2/J=l-0 ratio and CO J=3-2 intensity: 

R<1.7 0.568(0.024) -1.019(0.032) 0.89 

1.7<R<2.7 0.515(0.020) -0.878(0.019) 0.89 

2.7<R<3.7 0.549(0.019) -0.761(0.011) 0.80 

3.7<R<4.7 0.555(0.038) -0.683(0.016) 0.58 

4.7<R<5.7 0.551(0.022) -0.678(0.010) 0.68 

5.7<R<6.7 0.791(0.042) -0.630(0.013) 0.69 

R>6.7 1.104(0.104) -0.591(0.017) 0.46 
Relationship between PAH 8|Lim and CO J =3-2 intensity: 

R<1.7 0.869(0.044) 0.016(0.059) 0.85 

1.7<R<2.7 0.556(0.024) 0.285(0.019) 0.80 

2.7<R<3.7 0.643(0.017) 0.230(0.011) 0.91 

3.7<R<4.7 0.711(0.019) 0.202(0.010) 0.89 

4.7<R<5.7 0.805(0.021) 0.209(0.012) 0.88 

5.7<R<6.7 0.895(0.031) 0.249(0.013) 0.77 

R>6.7 0.989(0.049) 0.160(0.015) 0.45 

All 0.718(0.09) 0.220(0.005) 0.89 

Column 1: radial intervals (in kpc); Columns 2 and 3: parameters of linear 

fits to log(i^3_2/i_o) versus log(CO J=3-2 intensity) and log(PAH 8|j.m 

surface brightness) versus log(CO J=3-2 intensity); uncertainties are given 

in parentheses; Column 4: Pearson correlation coefficient. 



Table 7. Parameters of relationship in Equation]!] 



Radius 


A 


B 


C 


(kpc) 








R<1.7 


0.6 


0.19 


1.3 


1.7<R<2.7 


0.6 


0.12 


1.2 


2.7<R<3.7 


0.6 


0.12 


1.2 


3.7<R<4.7 


0.6 


0.12 


1.2 


4.7<R<5.7 


0.6 


0.12 


1.3 


5.7<R<6.7 


0.9 


0.10 


1.3 


R>6.7 


1.2 


0.09 


1.3 



Column 1: radial intervals (in kpc). The mean value in each interval is 
used; this is the central value in each radial interval, except for R< 1 .7 kpc 
where the mean radius is 1.1 kpc and R>6.7 kpc where the mean radius os 
7.6 kpc; Columns 2, 3 and 4: parameters of the linear fits to Equation[2j. 

brightness that was not apparent for the relationship with galacto- 
centric radius (FigurefTTJ. First, in the inter- arm regions, for a given 
CO line intensity the CO ratios R^-2/2-1 and i?2-i/i-o tend to 
be higher in Arm II than in Arm I. This difference is not apparent 
for the inter- arm in i?3_2/i-o, and is also not obvious for the spiral 
arms. Second, for R^_2/2-i and i?3_2/i-o the slope of the rela- 
tionship between CO ratio and CO surface brightness appears to be 
steeper for the inter- arm regions than for the spiral arms. Thirdly, 
for all three CO ratios we see that the relationship with CO sur- 
face density clearly varies with galactocentric radius, and that the 
dependence of the relationship on radius appears to be greater for 
R2-1/1-Q and i?3_2/i-o than for R^_2/2-i- At the centre of the 
galaxy, where the highest CO surface densities are located, there is 
little variation in the CO ratio for a given CO line intensity. How- 
ever, as galactocentric distance increases, and CO surface densities 
decrease, there is increasingly large variation in the CO ratio for a 
given /co- 
in Section |5] we noted an apparent trend for i?3_2/i-o in the 
inter-arm regions to increase with radius. From inspection of Fig- 
ure [12] (lower middle panel) we can now see clearer evidence of 



this. Here, for a given CO intensity high values of the ratio always 
occur at higher radii. This is in contrast to the behaviour seen in 
the spiral arms (lower left panel of Figure [T2|, where larger radii 
exhibit higher ratios for lower CO intensities than small radii, and 
likewise smaller radii exhibit low ratios for higher CO intensities 
than large radii. These trends can be explained as follows. The 
general dropoff in all CO ratios with radius, for both spiral arms 
and inter-arm regions, is likely explained by the decreasing beam- 
averaged molecular gas column density with radius. For the spiral 
arms, which contain star-forming regions, high ratio values reflect 
high gas densities and/or temperatures. The trend with radius sug- 
gests that active star-formation sites occur at lower CO intensities at 
larger radii, which may reflect an underlying star formation thresh- 
old (e.g. the Toomre Q parameter; |Toomre & Toomre|1972| ). One 
would not expect the same trend in the inter-arm regions since they 
contain little or no active star formation (see Figure|9] which shows 
that with the exception of inter-arm II to the north-east at large 
radii there are no significant regions of star formation that are not 
included within our spiral arm mask). There are a number of possi- 
ble explanations for the apparent trend with radius in the inter-arm 
regions. One possibility is that the temperature and/or density is in- 
creasing with radius, but this does not seem likely in the inter-arm 
regions. Another possibility is that the inter-arm medium consists 
of a diffuse molecular medium (dominating the lowest emission 
lines most) in which dense clouds are embedded. Loss of the dif- 
fuse medium at large radii would result in the observed behaviour. 
An alternative possibility is that the CO begins to become opti- 
cally thin at large radii. Under typical excitation conditions all low 
CO lines are optically thick, and if the density and temperature 
are high enough CO J =3 -2 will have higher optial depth than the 
lower lines. Thus, if cloud column densities in the inter-arm region 
decrease with radius, CO J=l-0 may become optically thin while 
CO J =3-2 remains optically thick, hence producing the observevd 
trend. 

Variation of the CO ratio with CO surface density suggests a 
relation to star formation and local heating. Variation with radius, 
on the other hand, suggests a relation to a more diffuse component 
in the disk. Thus, Figure [12] clearly suggests that the CO J=l-0 
line emission is also tracing a more diffuse and radius-dependent 
component of molecular gas than the higher CO transitions. 

While it is clear that the variation of CO ratio with CO surface 
density is dependent on galactocentric distance, the important result 
from this plot is that a simple relationship between the CO ratio 
i?3_2/i-o and the CO J =3-2 integrated intensity cannot be used 
to scale the CO J=l-0-derived X-factor for use in determining 
molecular gas surface density from the CO J=3-2 line without talk- 
ing into account the spatial location of the emission in the galaxy - 
the variation in the relationship implies that applying a global scal- 
ing for sub-kiloparsec scales will overestimate the molecular gas 
surface density at large galactocentric radii and underestimate it at 
small radii. 

In Figure [13] we investigate the radial dependence of the rela- 
tionship between CO ratio and CO J =3-2 intensity, by plotting the 
data at radial intervals of ^ 1 kpc and fitting a bi-variate linear func- 
tion for each bin (shown as solid lines in Figure [T3]l. The parame- 
ters of the linear fits are given in Table[6] We find that the Pearson's 
correlation coefficient decreases considerably with galactocentric 
radius, from 0.95 in the centre to 0.5 at ^7 kpc. 

The dashed lines in Figure [13] show radius dependent relations 
given by a linear function of the CO ratio with both CO J =3-2 
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Figure 13. Relationship between CO J=3-2/J=l-0 line intensity ratios and 
CO J =3-2 intensity at different galactocentric radii (over the entire galaxy 
disk where individual pixels are detected above 3cr in both CO transitions). 
We colour-code the plotted points according to galactocentric radius at 
smaller radial intervals than for Figure Q2] here we show the central emis- 
sion with R<1.7kpc and 1 kpc intervals thereafter. The solid lines show 
bi-variate least squares fits, whose parameters are listed in Table [6] The 
dashed fines show the radius-dependent relationship given by Equation [2] 
using the mean distance of pixels in each radial bin. 



intensity and galactocentric radius such that 



log (CO J = 3 - 2/J = 1 - 0) = A\og 



Ico 



Kkms 



+B 



B_ 
kpc 

(2) 
for the parameters given in TablejV] We find that a relationship with 
a slope of 0.6 reasonably describes the emission between radii of 
1.7 and 5.7 kpc, but that different relationships are needed to de- 
scribe both the central regions (shallower slope) and regions at radii 
larger than ^6 kpc (steeper slope). Nonetheless, a global fit across 
all radii yields a slope of 0.6, which is comparable to the value 
found by lBendoetaLlpOTOl for NGC 2403. The relationship be- 
tween the CO ratio and Equation |2] for individual pixels is shown 
in Figure [14] The scatter in the relationship clearly increases with 
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Figure 14. Relationship between CO J=3-2/J=l-0 intensity ratio and the 
relationship given by the right-hand side of EquationJ2] for all pixels de- 
tected above 3cr. Points are colour-coded as in Figure [13] 



radius, as indicated by the decreasing correlation coefficient (Ta- 
blelHl. 



6.6 Relations between CO J=3-2, PAH and 24|^m emission as 
a function of the spiral structure 

It is interesting to compare the relationships between PAHs, cold 
dust, molecular gas and atomic gas emission in M 51 with those in 
other late-type galaxies, and in particular whether PAH emission is 
a tracer of gas and dust in the cool ISM. A good target for a simple 
comparison is the nearby (^ 3 Mpc; |Freedman et al.|2001| ) SABcd 
spiral galaxy NGC 2403, whose relationships between these prop- 
erties were studied by |Bendo et al.|p010 ). While previous authors 
have found that PAH emission in nearby galaxies is strongly cor- 
related with cold dust emission (e.g. |Bendo et al.|2006| |Zhu et al!] 
|2008| ), and some authors (e.g.JRega^et al. 2006]) have found a cor- 
relation between PAH emission and CO line emission in nearby 
galaxies, Bendo etal. qOlO ) find that, in NGC 2403, the CO J=3-2 
and PAH surface brightnesses are uncorrected on sub-kpc scales. 

For both M51 and NGC 2403, as for many late-type spiral 
galaxies, the strongest regions of star formation are located outside 
the centre of the galaxy. Qualitatively, the structures seen in the 
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CO J=3-2, PAH, 2A\im and Ha images for M 51 and NGC2403 
are similar. Exceptions for M 5 1 are the inner part of Arm II, the 
south-west where CO emission is seen to trace the spiral arm but 
relatively little PAH emission is seen, and the very centre where 
there is a dip in the PAH emission but not CO (the latter differ- 
ence, however, may be due to the different spatial resolutions). 
Conversely, an exception for NGC 2403 is a hole in the CO emis- 
sion at the galaxy centre that is not seen in PAH emission (jBendol 
let al.|2010| ). 

6.6.1 Relationships between CO J =3-2 and PAH 8\xm emission 



Following [Bendo et al.] ( |2010| ) we used the CO J=3-2 integrated 
intensity and PAH8|j,m surface brightness images presented in 
Figures [T] and |3] to produce a map and radial profile of the 
CO J=3-2/PAH 8^Lm emission ratio (left-hand panel of Figure [15] 
and the top panel of Figure[T6] respectively). The two input images 
are in units of K km s~ ^ and MJy sr~ ^ , respectively. As before, only 
pixels detected above 3 a in the input images are shown. 

The CO J=3-2/PAH 8|^m ratio ranges from 0.04 to 3.5 across 
M 51, with a median value of 0.8 and standard deviation of 14.1. 

The CO J=3-2/PAH 8^Lm ratio appears to broadly follow the 
spiral structure, with the highest values in the spiral arms and the 
lowest values in the inter-arm regions. 

In the spiral arms we find a mean CO J=3-2/PAH 8|^m ratio 
of unity, and values of ^0.7 in the inter-arm regions. The ratio ap- 
pears to decrease with radius for the spiral arms, falling by a factor 
of about two between ^2 and 6 kpc. There appears to be no obvi- 
ous overall trend with radius for the inter-arm regions, but we note, 
however, that the presence of strong features in the inter-arm pro- 
file at ^A and ^6 kpc, which are similar to features in the spiral 
arms, suggests that the spiral arms traced in PAH 8|j,m emission 
may be broader than for CO emission. It is clear that the variation 
of the ratio of CO J=3-2 to PAH 8^Lm emission in M51 is very 
different to the profile that was found for NGC 2403 by B endo et 
[ar] ( |2010| ), who showed that in that galaxy the ratio increases with 
radius but decreases sharply in the central region; NGC 2403 typ- 
ically has lower ratio values than found for M 51. It seems likely 
that the CO J =3-2 and PAH S\im properties for M 51 are very dif- 
ferent from those for NGC 2403. While Bendo et aLl ( [20T0l ) showed 
that, for NGC 2403, the CO J=3-2 and PAH 8^m emission does 
not trace the same structures on small spatial scales, in the present 
work for M 51 we find that the CO J=3-2 and PAH 8|^m emission 
do trace the same structures on sub-kiloparsec scales, as is also 
clear from simple visual inspection of Figures [3}|5] 

The most striking feature is that the highest values of the 
CO J=3-2/PAH 8|j,m ratio are located ^1.5-2 kpc to the south- 
west, in Arm I (see Figure [15]). Here, ratio values are up to a fac- 
tor of two higher than typically seen elsewhere in the spiral arms. 
From inspection of Figure [3] it is evident that this region of en- 
hanced CO J=3-2-to-PAH 8|am emission corresponds to a region 
of reduced PAH 8|^m and 24|^m surface brightness. Here, there are 
no strong features apparent in the H a image, but an integral field 
spectroscopy study of nebular emission lines in the central regions 
by |Blanc et al.| ( |2009] l does show an H a emission feature at this lo- 
cation after extinction correction and subtraction of a contribution 
from diffuse ionised gas. This is also just south along the spiral 
arm from the location where we have noted the highest CO J =3-2 
and CO J=2-l arm/inter-arm contrast (Section[3]). Conversely, two 
obvious regions of low CO J=3-2/PAH 8|^m ratio occur at sym- 
metrically opposite locations ^6 kpc to the north-east (Arm I) and 
south-west (Arm II) that are coincident with the locations of bright 



peaks in the PAH 8|^m, 24|j,m and H a images that are not peaks in 
CO. 

In the top panel of Figure [Tt] we plot the relationship between 
CO J=3-2 and PAH 8|j,m surface brightness, colour-coded for dif- 
ferent radial intervals as in previous figures. We find a strong cor- 
relation, with a Pearson correlation coefficient of 0.89 and the best- 
fitting bi-variate least squares relation is described by 

log(PAH8|^m) = (0.718=b0.009)log(CO J=3-2) + (0.220±0.005) 

(3) 
. The largest scatter arises at the largest radii, both from the outer- 
most parts of the galaxy disk and from inter-arm region II at large 
radii. The extreme outliers correspond to pixels in the outer-most 
regions of the disk. The remainder of the scatter appears to arise 
from a more subtle dependence of the relationship on radius. We 
examine this further in Figure [T8]by fitting separate linear relation- 
ships to different radial intervals, whose parameters are given in 
Table [6] We find strong correlations (Pearson's correlation coeffi- 
cients of 0.8-0.9) in every radial interval with the exception of the 
largest radii (0.45 at R>6.7 kpc). We note that the extreme outliers 
in the R>6.7 kpc radial bin correspond to pixels in the outer-most 
regions of the disk. 

With the exception of the innermost bin, R < 1.7 kpc, the 
slope of the relationship increases with increasing radius. The steep 
innermost bin corresponds to the region of the nuclear disk which is 
rich in molecular gas and exhibits high CO J=3-2/PAH 8^Lm ratios 
over most of the region, as seen in the left-hand panel of Figure [T5] 
The outliers below the fitted relation in this bin (red points) corre- 
spond to the region of highest CO J=3-2/PAH 8^Lm ratios (and also 
the region of highest CO arm-interarm contrast) described above. 

The results in the present work support the idea that PAH 
emission is a good tracer of molecular gas, and that, for M51 
at least, the correlation between CO J=3-2 and PAH S\im emis- 
sion holds on sub-kpc scales and across different physical re- 
gions (whether in the nuclear region, spiral arms or inter- arm re- 
gions) in the galaxy. With the exception of one region of enhanced 
CO J=3-2/PAH 8|^m ratio, the dependence of the relationship ap- 
pears to be one of galactocentric radius. 

One explanation for the relationship between CO J =3-2 and 
PAH 8|j,m emission could be if the formation of molecular clouds 
occurs in regions with stellar potential wells l ILeroy et"aL]|2008| ). 
The CO emission then arises from the new molecular clouds and 
the increase in stellar emission in the potential wells heats the 
diffuse ISM, with enhanced PAH emission occurring in the re- 
gions where the strength of the interstellar radiation field is in- 
creased (e.g. [Bendo et al.|2 008^. If this is the case in M51, then 
the strong correlation between PAH and CO emission on sub- 
kiloparsec scales suggests that the two processes are occurring in 
physically similar locations. Further evidence of this is the similar 
behaviour of CO J=3-2 with 24|j,m emission, shown in the centre 
panels of Figures [T5] and \T6\ The relationship between CO J=3-2 
and 24|^m emission (lower panel of Figure [TT] again colour coded 
for different radial intervals) is described by the best-fitting bi- 
variate least squares relation 

log(24M,m) = (0.814 ± 0.010)log(CO J=3-2) + (0.257 =1= 0.006) 

(4) 
with a Pearson correlation coefficient of 0.89. The relationship ap- 
pears to vary with radius in a very similar way to that seen for 
PAH 8|j,m emission. Again, the extreme outliers in the R>6.7 kpc 
radial bin correspond to pixels in the outer-most regions of the disk. 
An explanation of the decrease of both the CO J=3-2/PAH 8|j,m 
and 24|j,m/PAH 8|j,m ratios with radius would be if the PAHs are 
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Figure 15. Ratio maps. Left: CO J=3-2 integrated intensity (Kkms~^) to R\H 8fxm surface brightness (MJy sr~^). Centre: CO J=3-2 integrated intensity 
(Kkms~^) to 24|Lim surface brightness (MJy sr~^). Right: PAH 8|j.m to 24 |j.m surface density (both in units of MJy sr~^). The two CO ratio maps are shown 
on the same colour scale. All maps were produced as described in Sections[2!2]and|6.6| 



excited not only by UV photons (from massive star-forming re- 
gions and hence traced by CO J=3-2 and 24|j,m emission) but also 
by softer optical photons which become relatively more important 
in the outer regions of the galaxy (see e.g. |Uchida et al.|1998[|Li &| 
|Draine|2002l ). 



6.6.2 The relationship between PAH 8\im and 24\im emission 

Previous authors have found that PAH emission does not trace star 
formation on sub-kiloparsec scales. We investigate this for the case 
of M51 via the relationship between PAH S\xm and 24|j,m surface 
density. It is thought that PAHs originate in the outflows of carbon- 
rich AGB stars and their descendants or in interstellar shocks (e.g. 
|Tielens|r2008| and references therein). 24|j,m emission traces re- 
cent dust enshrouded star formation, though may also include a 
component of diffuse emission that is not associated with recent 
star formation. We examine the relationship between PAH S\xm 
and 24|j,m emission in the right-hand panel of Figure [15] and the 
lower panel of Figure [16] which respectively show a map of the ra- 
tio of PAH S\im to 24|j,m surface density (using images that both 
have units of MJy sr~^) and the radial profile of this ratio. The 
median ratio is 0.78, with a standard deviation of 0.28, and val- 
ues for individual pixels range from <0.1 to 1.8. There are several 
things that are immediately striking. First, the surface density of 
PAH 8|am emission dominates over 24|j,m emission (i.e. PAH 8p,m 
to 24|j,m surface density ratios greater than unity) in the two outer 
spiral arms and in regions of prominent spurs bridging the two spi- 
ral arms. Elsewhere in the spiral arms there are are few pixels with 
values of about unity, but none within the central ^2 kpc. Second, 
apart from the two outer arms, regions dominant in PAH 8|j,m emis- 
sion appear to be concentrated to the south of the galaxy, away from 
the companion. Third, 24|j,m surface density is larger than the sur- 
face density of PAH 8^Lm emission at the locations of all the bright 
24 \im peaks, but there are nonetheless some 24 ^im peaks at whose 



locations PAH 8|am surface density appears to be dominate. It is 
also clear that the ratio of PAH 8^Lm emission to 24|j,m emission 
increases with radius (lower panel of Figure [T6|, which is contrary 
to the findings of |Bendo et al] pOTO] ) for NGC2403. Lastly, it is 
also striking that the variation of this ratio as a function of distance 
along the spiral arms and interarms is most similar to the radial 
profile of the H I emission (see Figure [To|. 



7 SUMMARY 

We have presented a new CO J=3-2 molecular line map of the 
nearby grand-design spiral galaxy M 5 1 obtained with HARP-B on 
the JCMT. This map, covering an area of ^14x 10 kpc and having a 
spatial resolution of ^600 pc, is the first CO map in the CO J=3-2 
transition to cover the full extent of the disk in M 5 1 . 

We have used the CO J=3-2 map together with CO J=2-l and 
CO J=l-0 data and other data from the literature to investigate the 
molecular gas properties of M 5 1 as a function of distance along 
the spiral arms, adopting a simple approach for defining the spiral 
structure. We have found the following results: 

• We detect CO J=3-2 emission over an area of ^21 x 14 kpc, 
across the entire disk of M51 out to its companion NGC 5195. 
The distribution of CO J=3-2 emission follows a clear two-armed 
spiral structure with broad (^1 kpc) arms, and is in close agree- 
ment with the distribution of CO J=2-l emission. We detect inter- 
arm CO J =3-2 emission everywhere within a central area of 
^ 14 X 10 kpc, as well as several spurs, or bridges, of CO emission 
connecting the spiral arms. 

• We find that the contamination of submm continuum by in- 
band CO J =3-2 emission ranges from up to 20% in the spiral arms 
to <10% in the inter-arm regions. 

• By projecting the CO integrated intensity maps of M 51 to var- 
ious edge-on orientations, we have shown that the major axis profile 
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Figure 16. Ratios as a function of distance along the spiral arms and inter-arm regions, starting from a radius of 1 .7 kpc. Top: CO J=3-2 intensity (in K km s~ ^ ) 
to PAH 8|xm surface brightness (in MJy sr~^ . Middle: CO J=3-2 intensity to 24 |Lim surface brightness. Bottom: PAH 8|j.m to 24 fxm surface density. Symbols 
are the same as for Figures [TQland[TT] 
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Figure 17. Relationships between CO J=3-2 intensity and PAH 8|j.m surface brightness (top) and CO J=3-2 intensity and 24fxm surface brightness (bottom) 
for different locations in M51: spiral arms, inter-arm regions, and the entire galaxy including the outer disk (from left), for all pixels detected above 3cr. 
Symbols are the same as for Figure [^ For emission over the entire galaxy (right-hand panels), both relationships have a Pearson correlation coefficent 
of 0.89. The soHd lines show a bi-variate least squares fit given by log (PAH Sp-m) = (0.718±0.009) log (CO J=3-2) + (0.220±0.005) and log (24 [xm) = 
(0.814±0.010) log (CO J=3-2) + (0.257±0.006). 



of a galaxy with bright molecular gas rich spiral arms can mimic a 
molecular ring for some viewing angles. 

• For the CO J=3-2 and CO J=2-l transitions there is a clear 
difference between the variation of arm and inter- arm emission 
with galactocentric radius, with the spiral arm emission decreas- 
ing with radius and the inter-arm emission relatively constant with 
radius. Conversely, for CO J=\-0 line and H I emission the arm and 
inter- arm regions appears to follow a similar trend with radius. 

• The contrast between spiral arm and inter-arm CO J=3-2 in- 
tensity decreases with galactocentric radius, falling by a factor of 
^-^3 over two kiloparsecs. There is a similar but less pronounced 
trend for CO J=2-l line emission, but for CO J=\-0 and H I emis- 
sion the arm-inter-arm contrast appears relatively constant with ra- 
dius. 

• At this spatial resolution, the CO J=3-2/J=2-l line ratio 
varies from ~0. 1-1.0, with average values of 0.5 over the whole 
disk, 0.5 in the spiral arms, and 0.6 in the central ^l^ In the inter- 
arm regions, values range from 0.2-0.6. 

• There is no clear evidence that the CO line ratios vary as a 
function of distance along the spiral arms, but there does appear to 



be a trend for CO line ratios in the inter-arm regions to increase 
with radius. 

• We find a strong relationship between the CO J =3-2 inte- 
grated intensity and the PAH S\im surface brightness. The ratio of 
CO J=3-2 to PAH 8|j,m emission decreases with radius, and we 
find a strong relationship between this ratio and the CO J =3-2 in- 
tensity that varies with galactocentric radius. 

• A similar but slightly steeper relationship is found be- 
tween CO J=3-2 and 24|j,m surface brightness, with the ratio of 
PAH 8|j,m to 24|j,m surface brightness increasing with radius. 
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Figure 18. Relationship between CO J =3-2 intensity and PAH 8|Lim sur- 
face brightness at different galactocentric radii (over the entire galaxy disk 
where individual pixels are detected above 3cr). We colour-code the plotted 
points according to galactocentric radius as for Figure [B] The solid lines 
show bi-variate least squares fits, whose parameters are listed in Table[6] 



REFERENCES 

Aalto, S., Hiittemeister, S., Scoville, N. Z., & Thaddeus, P. 1999, 

ApJ,522, 165 
Arimoto, N., Sofue, Y, & Tsujimoto, T, 1996, PASJ, 48, 275 
Bastian, N., Gieles, M., Efremov, Y N., & Lamers, H. J. G. L. M. 

2005, A&A, 443, 79 
Bell, T A., Roueff, E., Viti, S., & Williams, D. A. 2006, MNRAS, 

371, 1865 
Bendo, G. J., et al. 2006, ApJ, 652, 283 
Bendo, G. J., et al. 2008, MNRAS, 389, 629 
Bendo, G. J., et al. 2010, MNRAS, 402, 1409 
Blanc, G.A., Heiderman, A., Gebhardt K., Evans II, N.J., Adams, 

J. 2009, ApJ, 704, 842 
Buckle, J. v.. Hills, R. E., Smith, H., et al. 2009, MNRAS, 399, 

1026 
Calzetti, D., et al. 2005, ApJ, 633, 871 



Currie, M. J., Draper, P. W., Berry, D. S., Jenness, T., Cavanagh, 

B., & Economou, F. 2008, in Astronomical Data Analysis Soft- 
ware and Systems, ASP Conference Series, Vol. 394, eds. R. W., 

Argyle, P. S. Bunclark, & J. R. Lewis, 650 
Dumke, M., Nieten, C, Thuma, G., Wielebinski, R., & Walsh, W 

2001, A&A, 373, 853 
de Vaucouleurs G., de Vaucouleurs, A., Corwin, H. J., et al. 1991, 

Third Reference Catalogue of Bright Galaxies, Springer- Verlag, 

New York 
Egusa, R, et al., 2009, ApJ, 697, 1870 
Freedman W L., et al. 2001, ApJ, 553, 47 
Foyle, K., Rix, H. -W, Walter, R, & Leroy, A. 2010, ApJ, 725, 

534 
Garcia-Burillo, S., Guelin, M., & Cernicharo, J. 1993, A&A, 274, 

123 
Garcia-Burillo, S., Combes, R, & Gerin, M. 1993, A&A, 274, 148 
Gordon, K. D., Engelbracht, C.W, Rieke, G. H., Misselt, K. A., 

Smith, J. -D., T., Kennicutt, R. C, Jr. 2008, ApJ, 682, 336 
Guelin, M., Zylka, R., Mezger, RG., Haslam, C.G.T., & Kreysa, 

E. 1995, A&A, 298, L29 
Heifer, T. T. et al. 2003, ApJS, 145, 327 
Helou, G., et al. 2004, ApJS, 154, 253 
Hitschfeld, M., Kramer, C, Schuster, K.R, Garcia-Burillo, S., & 

Stutzki, J. 2009, A&A, 495, 795 
Ho, L., Pillipenko, A., & Sargent, W 1997,ApJ, 112, 315 
lono, D., et al. 2009, ApJ, 695, 1537 
Irwin, J. A., et al. 2011, MNRAS, 410, 1423 
Israel, R R, & Baas, R 2002, A&A, 383, 82 
Israel, R R, 2005, A&A, 438, 855 

Israel, RR, Tilanus, R.RJ., & Baas, R 2006, A&A, 445, 907 
Kennicutt Jr., R. C, Armus, L., Bendo, G. et al. 2003, PASP, 115, 

928 
Kennicutt JR., R.C. et al. 2007, AJ, 671, 348 
Kohno, K., Kawabe, R., Tosaki, T., & Okumura, S. 1996, ApJ, 

461, 29 
Koda, J., et al., 2012, ApJ, 761, 41 
Koda, J., et al., 2011, ApJS, 193, 19 
Koda, J., et al. 2009, ApJ, 700, L132 
Knapen, J. H., & Beckman, J. E. 1996, MNRAS, 283, 251 
Kramer, C, Mookerjea, B., Bayet ,E., et al. 2005, A&A, 441, 961 
Krist, J., 2002, Tiny Tim/SIRTIP USer's Guide, SSC, Pasadena 
Kuno, N. et al. 2007, PASJ, 59, 117 
Kuno, N. & Nakai, N. 1997, PASJ, 49, 279 
Kuno, N., Nakai, N., Handa, T, & Sofue, Y 1995, PASJ, 47, 745 
Leroy, A. K., Walter, R, Brinks, E., Bigiel, R, de Blok, W J. G., 

Madore, B., & Thornley, M. D. 2008, AJ, 136, 2782 
Li, A. & Draine, B. T. 2002, ApJ, 572, 232 
Lo, K. Y, Ball, R., Masson, T. G., Phillips, Scott, S., & Woody, 

D. R, 1987, ApJ, 317, L63 
Louie, M., Koda, J. & Egusa, R, et al., 2013, ApJ, 763, 94 
Matsushita, S., et al. 2004, ApJ, 616, 55 
Matsushita, S., Kohno, K., Vila-Vilaro, B., Tosaki, T. & Kawabe, 

R. 1998, ApJ, 495, 267 
Mauersberger, R., et al. 1999, A&A, 341, 256 
Meidt, S. E., Rand, R. J., Merrifield, M. R., Shetty, R., & Vogel, 

S. N. 2008, ApJ, 688, 224 
Meijerink, R., et al. 2005, A&A, 430, 427 
Mentuch Cooper, E., et al., 2012, ApJ, 755, 165 
Muraoka, K., et al. 2009, ApJ, 

Nakai, N., Kun, N., Handa, T., & Sofue, Y 1994, PASJ, 46, 527 
Nakai, N., Kuno, N. 1994, PASJ, 46, 527 
Nakai, N., Kuno, N. 1995, PASJ, 47, 761 



26 C. Vlahakis et al, 

Nikola, T., Geis, N., Herrmann, R, Madden, S. C, Poglitsch, A., 
Stacey, G. J., & Townes, C. H. 2001, ApJ, 561, 203 

Rand, RJ., & Kulkarni, S.R. 1990, ApJ, 349, L43 

Reach, W. T., et al. 2005, PASP, 117, 978 

Regan, M., et al. 2001, ApJ, 561, 218 

Regan, M., et al. 2006, ApJ, 652, 1112 

Sakamoto, K. et al. 1999, ApJ, 616, L59 

Schinnerer, E., WeiB, A., Aalto, S., & Scoville, N. Z. 2010, ApJ, 
719, 1588 

Schuster, K. R, Kramer, C., Hitschfeld, M., Garcia-Burillo, S., & 
Mookerjea, B. 2007, A&A, 461, 143 

Shetty, R., Vogel, S. N., Ostriker, E. C, & Teuben, R J. 2007, ApJ, 
665,1138 

Tielens, A. G. G. M., 2008, ARAA, 46, 289 

Tilanus, R. R J., & Allen, R. J. 1991, A&A, 244, 8 

Thornley, M. D., & Wilson, C. D. 1994, ApJ, 421, 458 

Toomre, A., & Toomre, J. 1972, ApJ, 178, 623 

Tully, R. B., 1974, ApJS, 27, 437 

Tully, R.B., 1988, Nearby Galaxies Catalog (Cambridge Univer- 
sity Press) 

Uchida, K. I., Sellgren, K. & Werner, M. 1998, ApJ, 493, L109 

Walter, R, Brinks, E., de Blok, W J. G., bigiel, R, Kennicutt, R. 
C, Jr., Thornley, M. D., & Leroy, A. K. 2008, AJ, 136, 2563 

Warren, B. E., et al. 2010, ApJ, 714, 571 

Wielebinski, R., Dunke, M., & Nieten, Ch. 1999, A&A, 347, 634 

Williams, et al. 1994, ApJ, 428, 693 

Wilson, C. D., et al. 2009, ApJ, 693, 1736 

Zhu, Y. -N., Wu, H., Cao C, Li H. -N. 2008, ApJ, 686, 155 

This paper has been typeset from a TgX/ MlgX file prepared by the 
author. 



